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Qj Abstract 

^' 

The results of an extensive study of B-meson properties in quenched lat- 

•i-H . 

^^ , tice QCD are presented. The studies are carried out in the static quark limit 

jrt ' where the b-quark is taken to be infinitely massive. Our computations rely on 

a multistate smearing method introduced previously, with smearing functions 
generated from a relativistic lattice quark model. Systematic errors arising 
from excited state contamination, finite volume effects, and the chiral extrapo- 
lation for the light quarks are estimated. We obtain continuum results for the 
mass splitting Mb^ — Mb^, = 86 ± 12{stat) lg(sj/si) MeV, the ratio of decay 
constants JbJIb^, = 1-22 ± 0.04(stoi) ± 0.02{syst). For the B-meson decay 
constant we separately exhibit the sizable uncertainties in the extrapolation 
to the continuum limit (a — > 0) and higher order perturbative matching. We 
obtain /s = 188 ± 23{stat) ± I5{syst) tf{extrap) ± U{pert) MeV. 



1 INTRODUCTION 

Heavy-light mesons are ideal systems for lattice QCD studies. For hadrons that contain a 
single heavy quark, the dynamics of QCD simplifies as mg — > oo. In this limit, the heavy 
quark is on shell and at rest relative to the hadron|l| and the QCD dynamics becomes 
independent of 77iQ|y|. Heavy quark mass dependence can be extracted analytically to 
produce an effective action for a static quark and the remaining light degrees of freedom g, |j] . 
In the continuum, the resulting heavy quark effective theory (HQET) makes transparent 
the symmetry S and scaling relations between systems which differ by heavy quark spin 
or flavor p, 0, Hi- Furthermore for flnite heavy quark masses, the effective action can be 
improved order by order in l/mq. The ttiq — > cx) limit may also provide other insights into 
QCD dynamics. In particular, heavy-light mesons are physical systems with a single light 
valence quark and fully relativistic QCD dynamics. It is likely that much can be learned 
about constituent quark ideas in this simple setting S. 

These theoretical developments have immediate physical applications for B physics. 
Since the b quark mass is significantly heavier than the other mass scales (the QCD scale 
and the light quarks masses) which enter into the dynamics of B hadrons, it is likely to be 
a good approximation to treat the b quark in the ttiq —> cx) limit within B hadrons. 

Recent developments in lattice gauge theory have led to the possibility of calculating 
the masses and decay constants of Bq {q=n, d, s) mesons from first principles (QCD) with 
enough accuracy to be of both phenomenological and theoretical interest. In particular, we 
consider the mass difference Mb^ — ^J^b^, the decay constant /b„ ^ and the ratio fB^/fs^- 
We will only consider the lattice action appropriate to the static limit {mi, — j oo) and hence 
all the results reported have corrections of order KqcD/'nib- A variety of other methods have 
been developed to study B mesons on the lattice. Methods for treating the b quark using 
nonrelativistic actionsB], a non-zero velocity formulation pO[ , and a generalization of the 
usual Wilson action which is not constrained to quark masses less than the inverse lattice 
spacing |11[ are being actively pursued. 

The present study encompasses a systematic analysis of data for Mb and /s at four 
different lattice spacings a (with associated beta values /3 = 5.7,5.9,6.1, and 6.3) and a 
variety of physical volumes (in lattice units 12'', 16"', and 20'^) for one fixed spacing {(3 — 5.9). 
For each case above at least four light quark mass values (kappa values corresponding to 



pion masses in the range 300 - 800 MeV) were studied. A uniform and consistent fitting 
scheme was used in all cases. This is particularly important for an accurate extrapolation 
to results at physical light quark masses (k = n^^d and k, — Kg), and for assessing the 
a-dependence of the results. 

In the static approximation, the heavy quark propagator is reduced to a straight timelike 
Wilson line, making it possible to calculate correlators of spatially smeared Qq operators 
without having to compute light quark propagators from smeared sources. Because of this 
simplifying feature, the heavy-light meson system is an ideal place to develop sophisticated 
operator smearing techniques. Such techniques are indispensable for the accurate calculation 
of fs and other B-meson properties. Until recently, most such calculations have relied on a 
more-or-less ad hoc choice of smearing functions (e.g. walls [p^, cubes [|3|, or exponentials 
0, or Gaussians ||l||). 

In the present study, we have applied a multistate smearing method introduced pre- 
viously in Ref. pq |. This analysis provides some significant improvements over previous 
investigations. First, we have made a serious effort to construct smearing functions which 
closely resemble the actual Coulomb gauge wave functions of the valence Qq system as mea- 
sured on the lattice. As reported in a previous paper p^ , the heavy-light wavefunctions 
from lattice QCD are reproduced with remarkable accuracy by a simple relativistic quark 
model (RQM) Hamiltonian which contains the static QCD potential extracted from Wilson 
lines in Coulomb gauge. In addition to being an interesting statement about QCD dynamics, 
the success of the RQM has a practical consequence which we will exploit here. The RQM 
Hamiltonian provides a simple and precise way of constructing orthonormal sets of realistic 
smearing functions for lattice heavy-light calculations. For each value of /3 and lattice size, 
the static Wilson potential is calculated from the gauge configurations and used in the RQM 
to generate heavy-light smearing functions. The only tuneable parameter in this procedure 
is the light-quark constituent mass parameter /i in the kinetic term of the RQM Hamilto- 
nian. In practice, this parameter was initially selected by measuring the lattice QCD ground 
state wavefunction and adjusting /i to give the best fit for the RQM ground state. In some 
cases, after an initial multistate fit to the heavy-light propagators, it was found useful to 
iterate the procedure with a more finely tuned value of the constituent quark mass /i, using 
the more accurate wave functions obtained from the multistate fit. 

In addition to this method for constructing smearing functions, another important inno- 



vation introduced in the present study is the fitting procedure used to extract information 
from the heavy-hght correlators. Starting with the wavefunctions from the first M S-wave 
states of the RQM, it is relatively easy to construct the entire M x M matrix of correlators 
among the corresponding smeared Qq operators, as well as the "smeared- local" correla- 
tors between each of these operators and the local Qq source. Along with the local-local 
correlator, these form an (M + 1) x (M + 1) matrix. This matrix contains far more infor- 
mation than just the smearcd-smeared and smeared-local correlators of any single smearing 
function. In particular, the matrix contains information about excited states, which, when 
properly exploited, allows an accurate extraction of ground state properties even at very 
short time separations, where excited state contributions are still large. The method we 
introduce to accomplish this employs a x^ minimization procedure to simultaneously fit the 
(M + 1) X (M+ 1) matrix of correlators to a sum of M exponential (pole) terms, representing 
the contribution of the M lowest lying heavy-light eigenstates. (In all fits, we exclude the 
local- local correlator, which, at short time separations, is not well fit by a few low lying 
states.) The matrix coefficient (residue) of each pole term factorizes and can be written 
in terms of an M + 1 component vector whose entries represent the vacuum-to-eigenstate 
matrix element of each smeared operator. In practice we have found an M = 2 fit to yield 
fairly accurate results for ground state properties. For all of the fits used, the x^ pcr degree 
of freedom was less than 1.3. The multiparameter fits were carried out using the CERNLIB 
minimization routine MINUIT. 

The improved control over systematic errors gained from the multistate fitting method 
allows us to better address a number of issues. In particular the dependence of heavy- light 
meson parameters on both the light quark mass and the lattice spacing are examined in 
detail. One of the difficulties with previous analyses which prevented accurate chiral and 
a ^ extrapolations was in the arbitrariness of the smearing procedure. It is clear that any 
ad hoc smearing function will have a substantial overlap with excited states. Typically one 
tries to optimize the smearing function (e.g. by adjusting the size of the cube) and to go far 
enough out in time that excited states have died away. The approximate equality of smeared- 
smeared and smeared-local effective masses, combined with some indication of an effective 
mass plateau, are the main criteria of success in this procedure. Unfortunately, the procedure 
is somewhat subjective and it is difficult to rule out large systematic errors due to excited 
state contamination. An attempt to reduce these errors by extracting results from larger 



time separations leads to a rapid deterioration of statistics. Moreover, because of the well- 
known signal-to-noise difficulty for the heavy-light propagator in the static approximation, 
[ p5| ITGJ , the problem of isolating the ground state becomes more difficult at smaller lattice 
spacing. As a result, extrapolation to a = is particularly problematic. Furthermore, 
use of any fixed smearing function at different light quark masses introduces a significant 
systematic error in the extracted kappa dependence. This is important in the determination 
of results for /bJIb^ and Mb^ — Mb^,- The multistate fitting procedure effectively deals 
with these difficulties, greatly reducing our errors. 

Extrapolating to the continuum limit (a = 0) we obtain the ratio Jb^/Jb^ = 1.22 ± 
.04(stat) ± .02(syst) and Mb, - Mb^ = 86 ± 12(stat) tglsyst) MeV for the mass difference. 
For these quantities, only a slight dependence on the lattice spacing is observed, and the 
systematic errors associated with the a ^ extrapolation (included in the above) are small. 

The situation for the decay constant is more complicated. We find a significant lattice 
spacing dependence for the ground-state pseudoscalar decay constant /s. The results for 
the four j3 values are consistent with either a linear or quadratic dependence on the lattice 
spacing a. The linearly extrapolated result in the a ^ limit is Jb = 188 ± 23 ± 15 MeV 
This result is notably smaller than previous estimates of /b in the static approximation. 
The primary reason for this is the a —t extrapolation. The quadratic extrapolated result 
in the a ^ limit is Jb — 214 ± 13 ± 17 MeV. This fit reflects the fact that our results at 
(3 = 6.3 are not inconsistent with those of Ref. |l^,[Q,||l^ (the latter two results using 
a clover action for the light quarks). In both cases, the first error is statistical while the 
second includes systematic errors, which we discuss in Section 5. Our final result for the Bu,d 
meson decay constant (in the heavy quark limit) /b = 188±23(stat)±15(syst) ]']Q^(extrap)± 
14(pert) MeV explicitly separates out this theoretical uncertainity associated with the a -^ 
extrapolation as well as our estimate of the uncertainity associated with higher order 
perturbative matching corrections. 

We will give a complete discussion of the perturbative matching of lattice to continuum 
results in Section 2. We include a discussion of the heavy quark mass renormalization in 
relation to the residual mass parameter (A = Mb — rnf, (quark)) of HQET. In Section 3 we 
discuss our analysis procedure. The details of the multistate smearing technique and the 
construction of the smearing functions from a relativistic quark model (RQM) are presented. 
In Section 4 the numerical lattice results at each /3 are presented. The statistical and fitting 



errors associated with our final physical results are dcterniincd. The light meson results used 
to set light quark masses are contained in Appendix A. The discussion of systematic errors 
associated with excited state contamination, finite volume, nonzero light quark masses, 
scale uncertainities, and the extrapolation to zero lattice spacing are all discussed in Section 
5. A study of the time evolution of the wavefunctions for heavy-light states is presented 
in Appendix B. These results provide an independent check that our multistate smearing 
analysis has removed excited state contamination. In the Section 6 we present our final 
results, compare them with other recent calculations, and discuss upcoming studies. 

2 Perturbative Matching 

2.1 Extracting Properties of Heavy-Light Mesons from LQCD 

In this section, we will focus on short-distance corrections to the results obtained from 
lattice QCD. These corrections are common to both traditional techniques for extracting 
meson properties and to the multistate smearing method employed here. Explanation of 
the details of the multistate smearing method are deferred to the next section. 

The corrections to the matrix elements of time component of the heavy-light axial 
current, J{)^{n), are computed by demanding that the ratio of the current renormalizcd 
with some continuum regularization scheme and the lattice-regularized current be unity. To 
be a little more precise, one computes this ratio, Z, using some matrix element, and from 
then on one multiplies any result obtained using the lattice-regularized current by Z . The 
states used to determine Z can be chosen for calculational convenience since the ratio is 
independent of the choice of states. Although there is no choice of states for which the 
numerator and denominator of the ratio are separately calculable, because the operators 
only differ at scales on the order of the cutoffs of the two regularizations, and at these scales 
QCD is perturbative, the ratio can be calculated in perturbation theory. 

The procedure is actually slightly more complicated than explained in the preceding 
paragraph. Because the lattice QCD calculations are done in the heavy quark effective 
theory (a theory which does not have the same particle content as the full standard model), 
it is necessary to compute an additional ratio, Z^ont- This is the ratio of the axial current 
renormalizcd at the scale m^ in the standard model to the axial current renormalizcd at a 
scale q* in the heavy quark effective theory. 



The calculations of Z and Zcont are the subject of the following two subsections. An 
analysis of the heavy quark mass renormalization follows thereafter. In the last subsection, 
we summarize the various constants used in this study. 

2.2 Calculation of Z 

The ratio Z introduced above was calculated some time ago, but there is substantial un- 
certainty in these calculations. Tadpole- improved perturbation theory, as formulated by 
Lepage and Mackenzie pO], promises to reduce these uncertainties below the ten per cent 
level at one-loop. The application of the tadpole-improvement program to heavy quark 
effective theory has recently been discussed by Bernard [pll, and calculations have been 
performed by Hernandez and Hill |22|| . In this subsection we summarize the calculation of Z 
within the framework of tadpole-improved perturbation theory. Hernandez and Hill consid- 
ered both the zero-separation and unit-separation point-split heavy- light axial currents. We 
will restrict our attention here to the zero-separation heavy-light current used in our Monte 
Carlo calculations. We will further restrict our attention to the case of Wilson fermions with 
r = 1 . The use of tadpole-improved perturbation theory results in a substantial reduction 
in our best estimate of the central value and uncertainty for ]b ■ 

The root of the tadpole-improvement program is a nonperturbative renormalization of 
the basic operators in the lattice action. These redefinitions absorb the large renormaliza- 
tions arising from lattice tadpole graphs. A related additional part of the Lepage-Mackenzie 
prescription is the use of a larger perturbative coupling. If one uses /3 to determine the per- 
turbative coupling, aiat, one-loop perturbative corrections are consistently underestimated. 
These perturbation theory problems are due to the fact that aiat is a poor choice of expan- 
sion parameter. For example at an inverse lattice spacing of 2 GeV, the tadpole-improved 
expansion parameter is ay = 0.16, which is twice as large as ajat- Lepage and Mackenzie 
argue that the best way to arrive at ay is from a non-perturbative lattice determination of 
a perturbatively calculable quantity, such as the gauge field plaquette expectation value. 

Using tadpole-improvement of the Wilson action for quarks on the lattice as a guide, 
one can perform tadpole-improvement of the heavy quark action, and this has been done in 
Ref. |2^ . Instead of discretizing 

S^ [ d'^xb^ido + gAo)b 



as 



S = ia^ J2„ ^^H ( ^(") - Uo{n-6ybin-6) J , 
it is discretized as 

Stadpoic-improvcd = »a^ X! ^^^'^^ ( ^("^ Uoin-6)'' b{n-6) 



Mo 

where uq is defined as 

The combination U^j,{x)/uq more closely corresponds to the continuum field (1 + 2agA^(a;)), 
than docs U^j,{x) itself. With the tadpole-improved action, there is an additional factor of 
1/uo for each gauge field link in the product. Thus the Green's function of two heavy-light 
currents separated by no lattice spacings in the time direction satisfies, 

\r ( \^ - ^BJno) 

[^ B\^Qj\tadpolc — improved — ^o 

The B meson decay constant /s is usually extracted from numerical simulations by 

fitting Gsino) to 

{fBrnsf r „ -, 

— exp[-CnoaJ 

2mB 

Thus we see that the tadpole improvement procedure has no effect on the fitted value of 
Jb- Its only effect is the change 

c^c + i^, 

a 
that is, a linearly divergent mass renormalization. 

So far we have seen that tadpole improvement does not affect the extraction of fs as 
it is generally done in lattice Monte Carlo calculations. However we must still take into 
account the effect of tadpole improvement of the light quark action, and this will involve 
some additional factors. 

As conventionally defined in lattice Monte Carlo calculations the lattice operator j" 
involved in calculating /^ is renormalized by a factor ^/2kcZ , where Kc is the critical value 
of the hopping parameter for the light quarks. The tree level value of Kc is 1/8. Lepage and 
Mackenzie advocate a reorganization of perturbation theory such that a factor of ^/Skc is 
included in Z and the renormalizing factor becomes Z /2. 

Let us see what this factor does at one-loop. Calculations of Z have been carried out to 
one-loop order and the result is of the form. 






Z=l + ^[ -4g{q) + ^Hq*a 



i2l 



This definition of g{q) (and similar ones for h{q),j(q), and k{q) which will be introduced 
below) follow the definitions in Ref. p^. A one- loop calculation of 8kc has been performed, 
and the result is expressible as 

1 = 1 ~ — /" —h( ) 

O'^c. onc—looD 'JTr ./ IT 



■^c, one— loop 

The relationship between Z and Z is 



Z — yoKc^onc—loopZ. 

The one-loop expression for Z is therefore 

5(g) = .9(9) + Hq)/^ 

where 

We continue with the application of the Lepage-Mackenzie prescription to determine the 
A- value of the coupling and the scale q* at which it is evaluated. The prescription for fixing 
the value of the coupling, ay is to extract it from a non-perturbative calculation of the 
1x1 Wilson loop {i.e., the expectation value of the plaquette, Upiaq). Once the coupling is 
known at some scale (alternatively, once the value Ay is known), it can be run to any other 
scale. The formula which relates ay to the non-perturbatively determined (lattice Monte 
Carlo) plaquette expectation value is, 

- ln(iTrC/p,,q) = y ay (3.41/a) [l - ay (3.41/a)(1.19 + 0.025n/) + 0{al)\ . 

The coefficient of Uf is the one appropriate for Wilson fermions with r = 1 . In the quenched 
approximation, Uf = 0. Ay is determined through 

av{q) = [/3oln(qVA^) + (/3i//3o) Inln(gVA'y))] "'. 

It remains to fix the scale q* at which av{q) is evaluated in the expression for Z. 
Lepage and Mackenzie propose to do that by calculating the expectation value of In q^ in 
the one-loop perturbative lattice correction. The formulae determining this scale are: 

Y = J d^q ~giq) ^ -13.93 (1) 

(ln(ga)2) ^ iy'A.9(<z)ln(gV) = |^ (2) 

q*a = exp[(ln(qa)2)/2] = 2.18 . (3) 



Hernandez and Hill quote the errors on the numerically evaluated values of Y and y In < 
q^ > as order 1 in the last decimal place. 

Using the two-loop formula for the running of ay with zero quark flavors one obtains 
av{q*)- The explicit dependence on the value of a^^ drops out of the ratio q* /Ay. Hence, 
the only way lattice Monte Carlo results have been used so far is for the expectation value 
of the plaquette; the determination of a~^ has not yet entered. The results for Z at various 
values of /3 are summarized in Table |l|. 

This completes our calculation of Z. It remains to multiply Z by the continuum running 
and matching factor, Zcont- We leave this for the following subsection. 

2.3 Continuum Running and Matching Factor 

For consistency, the one loop computation of Z should be combined with a two loop running 
in the continuum effective theory and a further one loop matching between the continuum 
effective and full theories. This produces the continuum running and matching factor Zcont 
which multiplies Z to give the full perturbative correction. It will turn out that the lattice to 
continuum matching factor Z is most significant, while Zcont produces only a small additional 
change. 

In the previous section, the coupling ay was determined in a no-flavor (quenched) lattice 
theory. We must now match onto a four-flavor (or five-flavor depending on the value of q*) 
continuum theory. In order for the infrared behavior of lattice and continuum theories to 
match exactly, we could use a continuum coupling whose value is equal to av{q*) (although 
differences between couplings are higher order effects). This implies we should choose a 
continuum scale qcont according to acontilcont) = c^v {<}*)■ In practice we ignore this criterion 
and simply set qcont — Q*, using ay(g*) everywhere in the matching. We then run in the 
continuum theory using four or five flavors depending on whether q* is greater or less than 
the 6-quark threshold mass. 

The 6-quark threshold itself is determined as follows. We assume the pole mass rribpoie 
is known and relate it to the MS running mass mi,{fi) according to p3| 

mtpoic = rribifi) [ 1 + -^^ - + InQi'^ /ml{fij) 

Setting fi — 771^ — jrLi,{m,l), we fix the threshold mass 777,^ by solving 



rubpoic = mb{ml) 
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Aacontjml) 



We use the usual two loop result for a, 



cont- 



47r / 6ilnln(^VA2) 



5oln(/xVA2) V feoln(/iVA2; 

with 

foo = 11 - 2n//3, h = 102 - 38n//3. 

Here n/ is the number of light flavors and we take A = 175 MeV for five hght flavors pj]. 
For other 7i/, A is fixed by demanding that acont be continuous. Applying this procedure 
with nib pole = 4.72 GeV p3 gives 



m* = 4.34 GeV. 



Now that q* and ml are fixed we use them in the combined two loop running plus one 
loop matching formula for Zcont ■ The one- loop anomalous dimension calculation |26[ for the 



heavy- light axial current was extended to two loops by Ji and Musolf |27 and the two loop 
result was confirmed by Broadhurst and Grozin |29|] . The same authors |3C|] also confirmed 
the one-loop matching calculation of Eichten and Hill p^ . The result for Zcont is: 



aoont^rn^y^'""^ 

acont {q*) ) 

. acont{rnl) - acontiq") 
An 



(4) 



70 (ji _ h 
2&0 \7o bo 



ciacontiml) 



47r 

In this equation 70 and 71 come from the anomalous dimension 7 of the heavy-light axial 
current in the effective theory, 

^cont , / ^cont \ 

7 = 70^^ +71 (^ ' 
while ci comes from the effective to full theory matching at one loop in the continuum, 



obtained by the method described at the beginning of section 2.1. This matching produces 
a contribution to Zcont of 

-, , a^cont 
1+Ci— . 

An 
Note that Ji and Musolf [ p7| quote ci with a sign error in the term which differs for vector 
and axial vector currents. The values of 7^ and ci are: 

70 = -4 

71 = -254/9 - 567r727 + 20n//9 

ci = -8/3 
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Table 1: Length scales and renormalization constants used in this paper. The values 
of a~^ for =5.7, 5.9, and 6.1 are taken from Ref. [43]. For /? = 6.3, a~^ is estimated 
from that at 6.1 by one- loop asymptotic freedom. Plaquette expectation values are 
taken from [20]. Z is the lattice to continuum renormalization factor for the axial 
current, Zcont is the factor relating the continuum heavy quark effective theory to 
full QCD, and Z^ is the overall renormalization factor used in previous discussions. 



(3 


a-i(GcV) 


(^Tr[/p,,q) 


Z 


^cont 


Za 


5.7 


1.15(8) 


0.549 


0.73 


1.00 


.63 


5.9 


1.78(9) 


0.582 


0.77 


0.96 


.65 


6.1 


2.43(15) 


0.605 


0.80 


0.94 


.68 


6.3 


3.08(18) 


0.623 


0.81 


0.93 


.68 



Using equation (^ we determine the Zcont values given in Table |[ These are then 
combined with the Z's obtained in the previous subsection to give the overall perturbative 
corrections Z listed in Table || Here, to conform with our previous notation ||2^ we also 
list the quantity 

Za - ^^ (5) 



2.4 Renormalization of the Heavy Quark Mass 

Before leaving the subject of short-distance perturbative correction we will compute the 
tadpole-improved estimate of heavy quark mass renormalization. 

In the continuum static limit neither the heavy meson mass (Mb) nor the renormalized 
heavy quark mass (rrib) appears explicitly, only the residual mass combination A = Ms — inb 
remains. The precise definition of the renormalized mass parameter in the dimensionally 
regularized heavy quark effective theory has been discussed by Falk, Neubert, and Luke 132]. 
Since the heavy quark mass renormalization is a linear divergence it vanishes in dimensional 
regularization with minimal subtraction. 

On the lattice, the heavy quark mass is renormalized in the static limit. This mass 
shift, 5m, is proportional to 1/a with a perturbatively calculable coefficient. Hence the 
mass, Meff, calculated for the ground state B meson in this lattice theory can be expressed 
as 

Meff = A-Sm; (6) 
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in terms the continuum residual mass A and the mass shift Sm. It is clear from Eq.(ph that 
the mass M^ff is linearly divergent as the lattice spacing a — > 0. However, it would appear 
that by measuring M^ff and removing the tadpole-improved one-loop mass counterterm, 
we have a determination of A. Hence we could obtain mh as defined in Ref. p2] from 
TUB- Unfortunately, Bigi et al |33| argue that non-perturbative effects ruin the preceding 
connection. Even if that is so, it is still possible to verify that the linearly divergent piece 
of Meff is correctly determined by the perturbative calculation of 5m. 

At one loop, the lattice heavy quark mass renormalization is of the form, 

Sm — — / — 5-A:(q) . 



The one-loop correction to the self-energy is linearly divergent and positive. The coefficient, 
(5™, of the mass counterterm is negative. The integrand k(q) is given by p4| , 

-| -1 

^sin2(q„/2) . 



fc(9) = \ 



As already noted the effect of tadpole improvement on the fitted value of the heavy 
quark mass is to change 

a 
This non-perturbative reduction of the mass is accompanied by a reduction in the countert- 
erm: 

Jm = --^- / -^k{q) . 

OTT a J TT^ 

where 

Hq) = k{q)+jiq) 

As in the calculation of the scale q* , we need to compute the expectation vlaue of In q^ in 
the integral of k(q). The formulae determining this scale are: 

X = f d'^qk{q) = 10.07 (7) 

(InM^) ^ ^/^'«^(«)M'Z'«') = ^ (8) 

p*a = exp[(ln((ja)2)/2] = 2.04 . (9) 

The calculational technique used is exactly that of Ref. |g^ , and the errors on the numerically 
evaluated values of X and X\n < q^ > are order 1 in the last decimal place. 
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Table 2: Tadpole improved mass counterterm for a static quark at the /3 values used 
in this paper. The associated gauge couplings av{p*) are also listed. The scale p* 
is 2.04/a (see text). 



(3 


avip*) 


aSm + In uq 


5.7 


0.228 


-0.394 


5.9 


0.191 


-0.339 


6.1 


0.170 


-0.307 


6.3 


0.156 


-0.285 



The values of av{p*) at the (3 values used in this study and aSm + Ihuq are presented 
in Table |. 

Two methods for the determination of A are possible. The first of these is to do the 
tadpole-improved perturbative subtraction just described. This results in a statistically 
independent result for A for each (3 value. The second method is to use the fact that the 
subtracted term is proportional to 1/a, and the physical value desired is independent of the 
lattice spacing. A two-parameter fit of the mass over the four values of /? has two degrees of 
freedom, and one of the parameters is the one we desire. The latter procedure ignores the 
running of the coefficient of the 1/a term. 

2.5 Comparison and Summary 

It is worthwhile at this point to compare the results in Tabic m with the widely used value 
of Za of 0.8, which does not benefit from tadpole improvement. Consider the results at 
(3 — 6.1. From the table, we find ZZcont/^/StKc to be 0.68. Consequently, tadpole-improved 
perturbation theory results in a reduction of the physical value of Jb by a factor of 0.68/0.8, 
i.e., a reduction of 18%. 

As an aside, we note that for the clover action — also termed "improved" , but in the sense 
that order a effects rather than tadpole graphs are being incorporated — the axial current 
has been renormalized ||35| , but tadpole improvement has not been applied to this operator. 

Now that perturbation theory has been reorganized to include tadpole corrections to all 
orders, we expect our one-loop calculation of the renormalization factor to be accurate to 
about 7%. This estimate of the magnitude of the two- loop corrections is obtained simply 
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by squaring the largest one-loop correction for the various values of r quoted in Ref. |23], 
for both discretizations of the axial current considered there. 

While the values for Zcont depend on the estimate of the lattice spacing, this dependence 
is weak; an increase by 10% in the estimate of the scale results in a reduction of at most 1.5% 
in the value of Zcont- Thus this source of uncertainty in the Z factor is negligible compared 
to both the estimated size of the two-loop corrections or the direct dependence of fs on the 
lattice spacing. 
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3 Analysis Procedure for Multistate Smearing 

3.1 Relativistic Quark Model for Heavy-light Systems 

The rapid deterioration |15l M of the signal to noise quahty of Euchdean correlators of 
heavy-light mesons at large Euclidean time makes the choice of an efficient smearing scheme 
essential if we wish to extract accurately the properties of low-lying heavy-light systems. 
In the multistate smearing approach previously introduced llq], the coupling of smeared 
bilocal Coulomb gauge operators to higher meson states was reduced by using smearing 
wavefunctions obtained from a relativistic quark model (RQM). The basic features of such 
a model are (a detailed examination of the connection of such a model with the full field 
theory in the case of the t'Hooft model can be found in |Q;see also,[p7|): 



(a) the use of a relativistic kinetic term \J'p^ + /i^ (with /i a constituent quark mass) for the 
kinetic piece of the Hamiltonian, and 

(b) a static confining potential Vir), which can be chosen to be the static interaction energy 
obtained from correlators of temporal Wilson lines in lattice QCD. 

The importance of relativistic kinematics in determining the shape of meson wavefunc- 
tions (with a light quark) was already implicit, if not clearly recognized, in the puzzling 
persistence of purely exponential falloff (~ exp— Cr) of hadronic wavefunctions, instead of 
the more rapid falloff one might naively expect in a confining model (~ exp —Cr^l'^ for a 
nonrelativistic particle in a linearly rising potential) . This exponential falloff is due to the 
nonlocal character of the kinetic part of the RQM Hamiltonian 

K{\r-T'\) = ^/p^TJi^exp^p■ {r-r)d^p (10) 

_/ |-9/2 g-pk-r-'l^ |r-r'|^C50 (11) 



r 



which implies that ^n(r) satisfying 

d^r'Kil r - / |)*„(r') + V^(r)*„(r) = £;„^„(r) (12) 

cannot fall exponentially faster than e^^*" (if it did, the integral over r' in ( |l2| ) would 
be dominated by r' ~ 0, giving an asymptotic behavior ~ e^^'' for the kinetic term, in 
contradiction with the assumed asymptotic behavior of {E„ — V{r))'^{r) for any V{r) with 
power growth). In other words, irrespective of the power rise of the confining potential, 
relativistic kinematics automatically smears out the wavefunction of a light quark over the 
Compton wavelength corresponding to the constituent quark mass. 
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In the static limit in which the heavy quark mass is taken to infinity, the relativistic 
Schrodinger equation (O) gives a single parameter fit (the constituent quark mass fi is 
the only adjustable parameter after the static potential has been measured on the lattice) 
to a complete set of orthogonal spin-independent wavefunctions corresponding to arbitrary 
radial and orbital excitations of the heavy-light system. To minimize lattice discretization 
and finite volume artifacts in the comparison of RQM and lattice Monte Carlo results, we 
have generated a set of lattice smearing functions by solving a discretized version of (|l2|), in 
each case on lattices of the same size as those used in the Monte Carlos, and in each case with 
the static potential determined from Wilson line correlators in the same gauge configurations 
used to extract our quenched QCD results. Namely, on each lattice and for each f3 value, we 
have extracted a full lattice static potential V{f) by measuring the correlator of two Wilson 
lines of time extent T (in Coulomb gauge-fixed configurations), and separated by a spatial 
distance r. The potential is then extracted by going out in Euclidean time T until the static 
energy £ — — ^In < W{0,T)W^f,T) >= V{f) stabilizes (for example, with (3 —5.9 on a 
16'^ lattice, this occurs for T > 5). The static potential extracted at various (3 values and 
lattice sizes is displayed in Fig |3.1| . 

The procedure used for generating lattice smearing functions from the RQM is as follows. 
We wish to obtain orthonormal lattice wavefunctions which are eigenstates of a lattice RQM 
Hamiltonian defined on a L^ lattice (with fjv' lattice sites): 

Hfp = Kpp + V{r)5pp> (13) 

p V ^ 

Such an eigenstate, in a channel of given orbital quantum number (S,P,D etc), will corre- 
spond to a pole of the resolvent applied to a source wavefunction ^^"^ of the same orbital 
symmetry: 

^ ^ II E(;^^)-'*^°H^")II - oo (15) 

r' 

For the starting source functions ^'^"^ one may take for example a monopole localized 
at the origin for S-states, a dipole for P-states, etc. After the energy E is tuned close to 
an eigenvalue En (until R is at least 10^ larger than the background value), a smearing 
eigenstate ^smcar(?^) is extracted by renormalizing the vector ^-h ^^^^ ^'^ ^™^ norm. The 
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inversion oi E — H is performed by the conjugate gradient algorithm, with the muhiphcation 
of the kinetic term done in momentum space using a fast Fourier transform. 

In most cases we have found it adequate to fit the RQM constituent mass by matching 
the IS wavefunction generated by the above procedure to the Coulomb gauge Bethe-Salpeter 
wavefunction obtained at a roughly fixed Euclidean time (corresponding to time slice 4 at 
/3=5.9). For the particular case of /3=5.9, k=0.159, on a IG'^ lattice, a more detailed fitting 
procedure was used to determine the optimal choice of quark mass jj, in order to fit the 

(15) 

meson wavefunction at various time slices. The mean square deviation of ^smcar from the 
measured quenched wavefunctions at various times T, for various ii, is displayed in Table y. 



Table 3: Mean square deviation of RQM smearing and LQCD heavy-light wavefunc- 
tions (xlO^), for (3 = 5.9, k = .159 as a function of time, /i is the constituent quark 
mass parameter in the RQM Hamiltonian. 



M 


T = 4 


T = 5 


T = 6 


T = 7 


T = 8 


0.05 


3.37 


1.99 


1.36 


5.65 


5.95 


0.10 


0.75 


0.42 


0.20 


1.86 


2.14 


0.12 


0.37 


0.45 


0.36 


1.02 


1.30 


0.15 


0.40 


1.06 


1.15 


0.39 


0.67 


0.20 


1.75 


3.30 


3.61 


0.72 


1.00 



From Table a we see that the optimal choice for the constituent quark mass varies in 
the range 0.10-0.15 if we fit to meson wavefunctions on time slices 4 to 8. We have chosen 
/i=0.12 as the best compromise for f3~5.9, k=0.159. With this parameter fixed, we have 
generated, by the procedure outlined following (E2|), IS, 2S, 3S and 4S smearing functions 
to be used in the multistate analysis described in the next section. For other /3, k values, we 
have usually used two smearing states only. The careful tuning of the quark mass performed 
here reduces to a very small level (4 percent or less) the coupling of the exact ground state to 
the higher smearing states, but will not turn out to be essential to the extraction of accurate 
masses and couplings for the ground state. A detailed discussion of the dependence of the 
results on the RQM mass parameter chosen for smearing is given in Section 5.1. 

As described in a recent article [|l^, the RQM gives a single-parameter fit to all the 
excited radial and orbital meson wavefunctions of our heavy-light system. After fixing /i 
by a match to the IS wavefunction, we have found remarkable agreement with measured 
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excited state wavefunctions (for example, the IP state, cf |T7|). This agreement suggests 
that this Ansatz accurately describes at least the valence quark sector of the full mesonic 
bound-state. 

3.2 Multistate Smearing 

Our object in this section is to outline a general procedure for extracting the maximum 
usable information from the multistate correlator matrix: 

C^Ht) ^ J2 *imU(^1 < I q{f, T)Q{0, T)Q{0, OW, 0) | > *Ll„(r-) (16) 

rr' 

where q{Q) are light (heavy) quark operators, and the ^smcar {a=l,2,...N) contain the set 
of orthonormal smearing functions obtained from the RQM as described in the preceding 
section. From a set of Nc decorrelated gauge configurations, we begin with a corresponding 
ensemble of Nc statistically independent C°'''{T) matrices, from which a standard deviation 
matrix a [T] can be obtained directly. In addition to the smearing wavefunctions of the 
relativistic potential model, the set {vI/('')} also includes the local source generating the 
desired heavy- light axial- vector matrix element for extracting Jb ■ Other types of smearing 
(cube, wall, etc) may also be included to facilitate an objective comparison with other recent 
calculations. In (|l^), the heavy and light quark propagators in each gauge configuration 
are computed in Coulomb gauge. As we are dealing with global color singlet states on each 
time slice (color sums are suppressed) C°^ is well-defined and non-zero. 
Defining states 

\¥\T>^Y. *ilar(^1Q(0, T)q{r. T)\Q> (17) 

r 

we have 

^^{t) = <$",T|$^0> (18) 

M 

= Yl e"^"^ < $" I " >< '^ I *^ > +0{e^^"+'^) (19) 

where the states | n > are exact eigenstates of the lattice Coulomb gauge transfer matrix 
with eigenvalues e~^". The remainder term of order q^^m+iT ^^y^ of course be small at large 
Euclidean time, but in addition should have a small prefactor to the extent that our smearing 
functions ^smcar(?^) (a=l,2,..M) do a good job in representing the valence quark structure 
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of the low- lying states, and to the extent that more complicated Fock states (containing 
extra quark pairs, real gluons, etc) are not too important. 

Next, define mixing coefficients (in our case, they are real): 

w"„ =< $" I n >=< n I $" > (20) 

Neglecting the exponential contamination of order q-^m+iT^ ^g ggg ^^g^^ -j-j^g multistate 
coupling matrix can be fit to an expression of the form 

M 

C"^''(r) = ^.;''„z;^e-^"^ (21) 

n=l 

Of course, we cannot hope to extract M independent time-dependencies with N < M 
smearing wavef unctions, so only N > M will be considered. Typically we shall extract 
the maximum information from the lattice data by picking N = M + 1 (the extra source 
function being the local current needed for the extraction of fs)- 
The fit is performed by a chi-square minimization of 

a,b T=T^ ^ ' 

with respect to the fitting parameters {w°„,-E„}, over a fitting range T< < T < T> in 
Euclidean time. The fit is performed on an ensemble of Nc jack-knife coupling matrices 
obtained by replacing each in turn of the Nc coupling matrices by the average matrix and 
reaveraging. We have chosen | $^ >= Jaxiai(O) | >, so the parameters u^j should 
be interpreted as lattice f-parameters for the ground and excited meson states. En as the 
corresponding masses, and v"'^ {o.— l,2,..,M) as mixing coefficients indicating the degree of 
overlap of the exact meson states with our RQM smeared states | $" >. The sum over a, b 
in ( P4 ) does not include the local-local correlator a = b = N, which is not well described 
by a sum over a few low-lying states. Note that this fitting procedure automatically gives 
the lattice f-parameters without the need to divide by the square-root of the smcarcd- 
smeared correlator as in the usual approach. Moreover, the ensemble of Nc parameter sets 
{w°„,i?„} obtained in this way can be subjected to a straightforward statistical analysis to 
determine the error in each of these parameters separately, correlations between parameters 
(e.g. between masses and f-parameters), and so on. 

Once the overlaps < $" | n > have been estimated by a best fit of C°'^{T), a smearing 
operator can be constructed which is guaranteed to contain at most one of the first N exact 
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meson states, thereby removing any other exponential time-dependence to the e^^"+i^ 
leveL Specifically, if eaia2..aN i^ the totally antisymmetric symbol in N dimensions, the 
smeared state 

|$^>=eaia...a„n^'''»l*"^>' A = 1,2,..M (23) 

is guaranteed (to the extent that we have accurately extracted the mixing coefficients 'u'^„) to 
contain only the exact meson state | A >, together with contaminations from the {N + l)'th 
excited state and higher. An effective mass plot of the usual kind can then be obtained for 
the A'th state by displaying (we use a smeared-local correlator to minimize noise) 

mUT)^\n-^^^ (24) 

where 

C'^{T)=<^^,T\¥^°''\0> (25) 

Of course, this plot will be most fiat for the ground state A = 1, where the relative expo- 
nential contamination is reduced to the level Q-iEN+i-Ei)T ^ ^^^ where small admixtures of 
lower-lying states cannot creep in to distort the effective mass plateau. The effective mass 
plots for the ground state ai (3 ~ 5.7, 5.9, 6.1, and 6.3 for various kappa values are shown in 
Figs.[ pT4 - t3.2[ . The solid line in each of these plots represents the ground state energy ex- 
tracted from the full multistate fit over the time window indicated by the length of the line. 
For each /3 the time window for the multistate fit was chosen to be over approximately the 
same interval in physical units, viz. about |/"i to 1/m. Noting that the splitting between 
the ground state and the second excited state is found, in our multistate calculations, to 
be around 800 to 900MeV, the choice of hfrn for the lower end of the time window should 
provide an exponential supression of excited states by at least a factor of 3. In addition, 
our careful tuning of the smearing functions should produce a relatively small coefficient 
for the higher excited states. The equality of the smeared-smeared and smeared-local effec- 
tive masses exhibited in the plots confirms our choice of fitting interval. A more complete 
discussion of systematic errors due to excited states is given in Section 5.1. 

The fitting formula (03) is easily generalized to allow a global fit to the data at various 
K values (for fixed f3) : this is essential in order to take into account correlations between the 
coupling matrices at different k values, which would affect our estimate for the error of the 
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results when linearly extrapolated to Kc- Namely, ( p2| ) is replaced by 

a,b K T=T. 

where 



x-^EEE ' "'""^"irj";-"" (26) 



a''''(K,Ty 



v'-'^in) 


= 


V^On+^lni^-'- 


- -c') 


En{n) 


EE 


Eon+Ei„iK-^~ 


'-c') 



and the mixing coefficients w"^^'"*^ are varied freely. Note that only the immediately 
physical mass and lattice-f parameters are assumed to have the chiral dependence on k"^: 
other mixing coefficients involve the model dependent choice of smearing functions from the 
RQM. The chi-square minimization allows the direct extraction of masses and couplings 
extrapolated to Kc, as well as the slopes in k^^ of these quantities (all of which are free 
variational parameters in this new global fitting procedure). The usual jackknife procedure 
can then be applied to yield the correct errors on the extrapolated quantities. 
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4 Lattice Results for Heavy-Light Mesons 

To extract results for masses and decay constants we have used the set of gauge config- 
urations and Ught quark propagators enumerated in Table |j. The light quark action we 
use is not 0{a) improved. The four columns in Table ^ are the gauge coupling, /?; lattice 
size; number of gauge configurations (separated by 1000, 2000, 4000, and 4000 sweeps for 
(3 — 5.7, 5.9, 6.1, and 6.3 respectively); and the light quark k values calculated for each 
configuration. 

Table 4: Summary of gauge configurations and light quark parameters used in this 
paper. Listed are a letter used to identify each Monte Carlo run, the /3 value, lattice 
size, number of gauge configurations in each ensemble, and the values of light quark 
hopping parameter n analyzed. 



run 


P 


lattice 


confs 


K 


b 


5.7 


123 X 24 


100 


.168,. 1667,. 165, .161 


e 


5.9 


12-'^ X 24 


100 


.159, .158, .157, .154 


c 


5.9 


16^ X 32 


100 


.159, .158, .157, .156, .154 


f 


5.9 


20^ X 40 


100 


.159, .158, .157, .154 


d 


6.1 


243 X 48 


50 


.1545,. 154, .153,. 151 


g 


6.3 


323 X 48 


50 


.1515,. 1513,. 1510,. 1500 



The multistate smearing analysis outlined in Section 3 provides a powerful method for 
extracting heavy-light meson parameters. Unlike single or double exponential fits to single- 
channel "smeared-smeared" and "smeared-local" correlators, the fitting of the Nx N matrix 
of correlators to an expression of the form ( pi| ) is highly constrained. As we will show, this 
method allows a determination of fs and other heavy-light parameters which is less prone to 
systematic errors than previously applied methods. In this section, we present our results. 
In Section 6.2 we compare these results with those recently reported in Rcfs.]!^, |lj, |9| |41[ |. 
Most of the results presented here were obtained from the fitting procedure discussed in 
Section 4 using iV = 3 and M — 2, i.e. a 3 x 3 matrix of correlators (2 smearing functions 
and the (5-function source) fit to the sum of two exponentials (always excluding the local- 
local correlator from the fit). To estimate the systematic errors associated with the fitting 
procedure, we have tried varying both the shape of the smearing functions (by adjusting 
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the RQM quark mass parameter) and the number Af of smearing functions included. These 
results are discussed below. 

First consider the mass eigenvalue Ei in Eq. (|2^) which describes the leading exponential 
falloff of the heavy-light correlators 

C'^'^it) ~ <wje-^i* (27) 

In the multistate fitting procedure, Ei is the energy associated with the ground state contri- 
bution to the correlator. This parameter represents a combination of the binding energy of 
the B-meson ground state plus a divergent mass shift of the heavy quark. Recall that only 
the bare mass of the heavy quark is removed in reducing to the effective static theory. The 
mass shift induced by QCD is therefore measureable on the lattice as a non-scaling piece 
in the parameter Ei . We measure the dependence of Ei on both the light-quark hopping 
parameter k and the lattice spacing a. The graph in Fig. shows the ^-dependence of Ei 
for the four values of /3 studied. In each of the four data sets, the dependence on k~^ is 
quite linear, allowing an accurate extrapolation to the chiral limit Kc. Numerical results 
are tabulated in Table M For each value of beta and kappa the results in this table were 
obtained from a 2-state fit, with the x^ per degree of freedom of each fit listed in the last 
column. For /3 — 5.9, the x^/dof are those obtained on the 20'^ box. The results for /b 
and aEi at /3 = 5.9 listed in Table g are the infinite volume values obtained by fitting all 
three box sizes to a Luscher finite volume formula, as described in Section 5.2 (except for 
K — .156, which has only been done on the 16'^ box). The results for the three box sizes 
are listed separately in Tables and ^ A measure of the overlap between the true ground 
state and the RQM wavefunction smeared operators is given by \/2k X]a=i(''^i)^ /V6- 
The measured values of this overlap for each of the fits is recorded in Table ||. 

In addition to analyzing the data at each value of k separately, we have also performed a 
simultaneous fit to all n values for a given f3 by allowing the multistate fitting parameters u° 
to depend on k, as described in Section 3. Using this procedure on jackknifed subensembles 
provides a better estimate of the error on the chiral extrapolation, since it takes account 
of the fact that the different k values have correlated errors. Comparing the results of this 
analysis with results of a separate analysis of each k, we find that the extrapolated values 
for El and /b change very little, while the error bars on these results are about 30% lower 
than those obtained by ignoring n correlations. On the other hand, the results for the 
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Table 5: Lattice results for heavy-light mesons (static approximation). Values for 
the ground state energy aEi and decay constant fs are extracted from a 2-state 
fit over time window AT. (fs is related to the physical decay constant fs by Eq. 
(30)). Results at K = Kc are from the multi-K fits as discussed in Sec. 3.2. The 
column labeled overlap is a measure of the total overlap between the true ground 
state and the RQM wavefunction smeared operators used. 



/3 


AT 


K 


aEi 


overlap 


Ib 


xVdof 


5.7 


2-8 


.161 


.827(6) 


.725(7) 


.670(19) 


.50 






.165 


.794(8) 


.717(8) 


.626(23) 


.78 






.1667 


.776(9) 


.699(7) 


.590(24) 


.54 






.168 


.767(11) 


.694(9) 


.578(29) 


.43 






Kc 


.758(10) 


.694(8) 


.564(28) 


.58 


5.9 


3-10 


.154 


.719(5) 


.742(12) 


.347(11) 


.87 






.156 


.692(8) 


.731(13) 


.318(13) 


1.06 






.157 


.678(6) 


.724(13) 


.300(11) 


.62 






.158 


.665(7) 


.716(15) 


.283(12) 


.56 






.159 


.645(9) 


.704(18) 


.259(14) 


.60 






Kc 


.638(9) 


.686(24) 


.250(14) 


.66 


6.1 


4-12 


.151 


.620(7) 


.769(17) 


.199(10) 


.57 






.153 


.583(9) 


.721(24) 


.170(12) 


.60 






.154 


.561(11) 


.705(24) 


.149(12) 


.61 






.1545 


.551(13) 


.700(30) 


.142(14) 


.48 






Kc 


.544(12) 


.689(25) 


.135(13) 


.55 


6.3 


5-14 


.1500 


.528(7) 


.748(18) 


.120(7) 


.72 






.1510 


.511(7) 


.728(20) 


.107(7) 


.66 






.1513 


.506(8) 


.724(20) 


.104(7) 


.61 






.1515 


.504(8) 


.729(20) 


.103(8) 


.56 






Kc 


.499(9) 


.720(17) 


.099(8) 


.62 
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slope of these quantities as a function of k~^ are greatly improved by the simultaneous-K 
fit, reducing the errors by a factor of 3 or more over the independent-K analysis. Thus it is 
especially important to take account of inter-K correlations for quantities such as ttib, —friB^ 
and JsJfB^- 

Let us consider Ei as a function of the lattice spacing a and of the naive light quark 
mass (see Appendix). The linear dependence on k^^ becomes 

Ei{2amq, a) = Ei{Q, a) + (k^^ - k-^)E[{0, a) (28) 

= Ei{Q,a) + 2amqE[{0,a) 

The quantity M^fj defined in Section 2, Eq. (||) is just Ei evaluated at k = Kc, 

M,ff^EiiO,a) (29) 

The slope parameters E[ (0, a) are obtained by a correlated fit to all k- values, as described 
in Section 3. The results for the slopes and intercepts are given in Table o. 

Table 6: Slopes and intercepts for aEi and /b as a function of k~^ — k~^. For a 
given /3, values are obtained from a simultaneous two-state fit to all values of k over 
a time window AT. 



p 


AT 


aTi(O) 


Em 


MO) 


/s(0) 


5.7 


2-8 


.758(11) 


.236(25) 


.564(29) 


.424(72) 


5.9 


3-10 


.638(9) 


.350(25) 


.250(13) 


.444(38) 


6.1 


4-12 


.544(12) 


.450(47) 


.138(13) 


.387(46) 


6.3 


5-14 


.499(9) 


.376(56) 


.099(7) 


.290(56) 



In a similar way, we obtain a linear fit (see Fig. |j) to the n^^ dependence of the 
groundstate pseudoscalar decay constant fs- Define a quantity fs which is just the matrix 
element parameter v^ in the multistate fit (overlap of the ground state with the (5-function 
source), with a normalization factor V2k for the light quark included. In the scaling limit, 
the physical value of /s is related to /b by the following multiplicative constants. 



Ib = Ib-x 



Mf 



X a 



xZa 



(30) 
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where Za is the renormaUzation factor associated with matching the fuU theory with the 



effective static theory on the lattice, as discussed in p3, |42 , and in Section 2. 

The lattice spacing dependence of E\ in the chiral limit m^ = is plotted in Fig. 0. 
The results are consistent with a linear a dependence, 

a£;i(0,a) = £;i(0,0) + i^ia (31) 



with 

£;i(0,0) = .351(14) (32) 

Ex = .481(25)GeV (33) 

The first term £'i(0, 0) is the linearly divergent (i.e. 0(l/a)) term in the heavy quark mass 
shift. 

The tadpole-improved estimate of abm was discussed in Section 2. The corresponding 
quantity is a(5Tn + lnuo, and was tabulated in Table 0. Comparing the ground state effective 
mass with the tadpole improved 1-loop result at each lattice spacing one finds about a 
30 % discrepancy in the singular part of the mass shift, which can easily be accounted 
for by higher-loop and/or nonperturbative contributions. In fact, ordinary (non-improved) 



one-loop perturbation theory |34 gives 






X -^ X 19.95 (34) 



Simply identifying this value with the extrapolated lattice result gives as — 0.162(6) which 
is in reasonable agreement with other determinations of ols in the range of lattice spacings 
considered here Ig^l • So the entire "discrepancy" can be removed by a reasonable redefinition 
of the perturbative coupling being used. 

The mass of the meson i?s, composed of a 6 quark and a strange antiquark, is of 
considerable phenomenological interest. Our calculation of the heavy-light ground state 
energy as a function of k provides a determination of the mass splitting between the Bs and 
the i3„ mesons 

AMs^ = Mb, - Mb^ (35) 
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Since the divergent self-mass of the heavy quark is independent of hght quark mass, it will 
cancel in the mass difference AMb^ , and the latter should therefore scale properly with a. 
For each value of /?, we use the determination of Kg and «„ discussed in the Appendix, along 
with the observed k dependence of Ei to determine AMb^ ■ The results are shown in Fig. 



A linear extrapolation of the mass difference AMb^ to a = gives 

AMb, = 86 ± 12 MeV (36) 

Notice that the results for AMb, shown in Fig. ^ exhibit a fairly mild dependence on the 
lattice spacing, in marked contrast to the strong a-dependence of /s. The decay constant 
/b^ for the strange B-meson may be determined, using the values for the slope parameter 
f'g in Table g. The ratio /bJ/b^, is plotted in Fig. H. Again, the a-dependence of the ratio 
is much weaker than that of Jb itself. Extrapolating to a = 0, we obtain 

^ = 1.216 ±0.041 (37) 



28 



5 Determination of Systematic Errors 

5.1 Systematic error due to excited state contamination 

The main results presented for /b have been obtained by a 2-state fit to the correlators 
which employed quark-model smearing functions for the ground state (f S) and first radially 
excited (2S) state. As pointed out in Section 3, this produced quite stable effective mass 
plots which indicated that accurate ground state parameters could be extracted with time 
separations as short as T = 2 or 3. Since the errors in most previous calculations have 
been dominated by the systematic effect of higher state contamination, it is particularly 
important to estimate the size of this effect to get an overall determination of the accuracy 
of our results. To further investigate this issue, we have carried out a more complete study 
of the dependence of the extracted fs value on the fitting procedure. First, we have varied 
the size of the source smearing function by changing the quark mass parameter /i in the 
RQM wavefunctions. We then compare the results from the 2-state fit with those from a 
truncated f-state fit (using only the correlators of the ground state smearing function and 
the delta- function source). 



In Figure |3.l| we compare the effective mass from the I-state fit with that of the 2- 
statc fit. After determining the IS-2S splitting A from a 2-state fit, the effective masses 
obtained over a Euclidean time window can be plotted versus the variable e^ , allowing 
an extrapolation to t = oo (See below). The results of the I-state fit are plotted for 
four different time windows, 1-6, 2-7, 3-8, and 4-9, and for four different choices of RQM 
smearing functions (with quark mass /i = .32, .60, .90, and 1.20). (Note: The result from 
window t< — t> is plotted at the value of e^'^* corresponding to t = i<.) The effective 
masses from the 2-statc fit using the four different fj, values and the time window 3-IG are 
all plotted on the far left side of the plot at e^'^* = 0. The I-state results for the different 
smearing function choices are clearly converging to a common effective mass at i = oo 
which agrees well with the 2-state result, the latter being quite insensitive to the choice of 
/i. Similarly, the result for fg from the full 2-state fit remains unchanged, within errors, 
for a wide variation of the fi parameter. On the other hand, the result from the I-state fit 
varies by 20-30 % over the same range of /i values. This provides strong evidence that the 
2-state fit does a good job of isolating the ground state, even when the chosen smearing 
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functions are not very well tuned. To demonstrate this, we look at the dependence of these 
results on the time window chosen for fitting. While the two-state fits are generally stable 
under variation of the time window, the 1-state fit shows a systematic time-dependence. If 
this time-dependence is assumed to be largely due to contamination from the first excited 
state, it should fit asymptotically to the functional form 

fB{t)^fB{^) + Ce-^' (38) 

where A = £'2 ^ -^i is the energy splitting between the ground state and first excited 
state. Without an independent estimate of this splitting, it is difhcult to obtain a reliable 
determination of the parameters in ( pq ) directly from the results of the 1-state fit. On the 
other hand, the 2-state fit determines both Ei and £'2, and therefore A. If we use this 
determination to fix A, the formula ( p8| ) can be used to extrapolate the 1-state results to 
t = 00. Comparing this result with that of the full 2-state fit provides a useful and nontrivial 
check on the assertion that the systematic effect from excited states is under control. Figure 



5.1 shows the results of such a comparison for the case /3 — 5.7, n — .161. The splitting 



obtained from the 2-state fit is A = .321 in lattice units. It is seen in Fig. 5.1 that, for 
the time window 1-6 (far right on the graph), the result for fs varies systematically with 
the choice of smearing function. As t gets larger, the results from the different smearing 
functions tend to converge to the same value. The points plotted at e^'^* = include the 
four extrapolated values obtained from Eq. (|3g). For comparison, the results of the 2-state 
fit for the four /i values and i = 3 — 10 are also plotted. All of these points are well within 
a standard deviation of each other. Similar comparison of the results of 1-state and 2-state 
fits for other values of /? and k give comparable agreement. Based on this agreement, we 
conclude that the systematic error on our results due to excited states has been eliminated 
at the level of our present statistics. 

It is worth emphasizing here that our ability to control excited state contamination 
depends crucially on the use of the multi-state fitting procedure. Although the 1-state 
fits were all found to lead to consistent results after extrapolation to t = cx), an accurate 
extrapolation would not have been possible without an independent determination of the 
splitting A, which is only obtainable from the 2-state fit. 
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5.2 Finite volume corrections 

Using the scales in Table || we find that the physical volumes of the boxes for the main 
ensembles used in our calculation are approximately {2.0fm)^, (1.8/m)'^, (1.9/m)'^, and 
(2.0/m)3 for (3 = 5.7(123), 59(163)^ 6.1(243), ^^^ 6.3(323) respectively. Although these 
volumes appear to be comfortably large compared to the observed size of the ground state 
B-meson on the lattice, we consider in this section the possibility of corrections to our results 
due to finite volume effects and describe our method for estimating these effects. Although 
our overall conclusion is that these effects are negligible on the lattices considered here, 
the estimates discussed in this section may be useful for selecting /3's and lattice sizes in 
subsequent studies. 

An extensive theoretical study of finite volume effects on field theoretic calculations has 
been carried out by Luscher |4^. Consider, for example, the effect on the mass of a particle 
mp. For large enough volume, the leading effect is due to the propagation of the lightest 
mass meson (e.g. pion) "around the world," leading to the expression for m{L), the particle 
mass in an i X L X i box, 

g-AL 

mp{L) = mp{oo) + A—— (39) 

where the exponent A is determined by the mass of the pion, and A is given in terms of 
the on-shell nPP coupling. This finite volume correction can be interpreted as the effect 
of squeezing the pion cloud surrounding the particle. A somewhat different situation takes 
place when the particle P is a loosely bound state of constituents. In this case, the finite 
size effect is caused by the squeezing of the bound state wave function |M . As pointed out 
by Luscher W7|, this situation falls into the same general framework as that which led to 



Eq.(39), except that, in this case, the particle that travels around the world is one of the 
constituents of the bound state. In fact, for a nonrclativistic bound state in a non-confining 
potential, the finite volume effect assumes exactly the same form as (p9|), but in this case, 
the exponent A is related to the binding energy (and hence to the spatial extent of the bound 
state wave function). 

For the case (3 = 5.9 we have carried out a complete Monte Carlo investigation of the 
heavy-light propagators on lattices of three different sizes, 12^ x 24, 16^ x 32, and 20^ x 40 
(runs e, c, and fin TableQ). With the value a^^ — 1.78 GeV, these three boxes are of spatial 
length 1.3 fm, 1.8 fm, and 2.2 fm, respectively. The results for both the ground state energy 



31 



aEi and for /b are given in Tables M and pi They are seen to be the same, within errors, 
in aU three size boxes, and thus, no significant finite size effect is observed. In order to 
determine an upper hmit on the finite volume corrections to our results, we will make the 
assumption that these effects can be parametrized in the Luscher form (p9). 

„-XL 



aEi (L) = aEi (cx)) + AaEi ■ 



fB{L) = fB{^)+A. 



L 

XL 



(40) 



(41) 



/« L 

(Note: For particle masses in full QCD in a sufficiently large box, such an expression has 
been derived rigorously. For masses in quenched approximation and for decay constants, 
it's validity is not established, but we adopt it as a convenient ansatz. An alternative power 
law form is also discussed at the end of this section.) For a given choice of the exponential 
parameter A, a fit to Eq. (p9[) gives a limit on the coefficient A. 



Table 7: Volume dependence of aEi at /3 = 5.9. Results are obtained from a 2-state 
fit over time window AT = 3 - 10 on lattices of size 12^ x 24, 16^ x 32, and 20^ x 40. 
Numbers in square brackets are the x^ P^r degree of freedom for each multistate fit. 



K 


aEi{12) 


a£;i(16) 


a£:i(20) 


.154 


.716(9) 


.718(7) 


.719(6) 




[.74] 


[1.12] 


[.87] 


.157 


.675(12) 


.679(8) 


.677(8) 




[.92] 


[1.01] 


[.62] 


.158 


.661(14) 


.667(9) 


.662(9) 




[.99] 


[.92] 


[.56] 


.159 


.652(21) 


.652(12) 


.641(11) 




[.98] 


[.76] 


[.60] 


Kc 


.638(17) 


.643(11) 


.634(11) 




[1.01] 


[.92] 


[.66] 



Our strategy is to extract an estimate of the exponent A in the Luscher formula by 
two methods: (1) a direct study of the finite volume effects in the relativistic quark model, 
and (2) a study of the exponential falloff of the ground state wave function obtained in 
the LQCD calculation. For both f = Jb and / = ai?i, the results from the RQM were 
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Table 8: Volume dependence of /^ at /? = 5.9. (See caption of Table 7.) 



K 


fB {12) 


hm 


/b(20) 


.154 


.341(21) 


.347(13) 


.346(14) 


.157 


.299(22) 


.303(14) 


.298(15) 


.158 


.284(24) 


.288(15) 


.278(16) 


.159 


.279(35) 


.271(18) 


.252(17) 


nc 


.261(29) 


.260(18) 


.245(18) 



calculated on 12^,16'^, and 20'^ boxes. (The relativistic Van Royen-Weisskopf formula ||37| 
was used to obtain Jb from the RQM wave function.) For the lightest quark mass studied, 
the RQM estimate gives A/a = 0.9 GeV. (The results for aEi and for /b are both well fit 
with the same value of A.) A slight increase in the value of A for larger light quark mass is 
observed, but is inconsequential for our analysis. A direct study of the exponential falloff of 
the LQCD ground state wave function (using the bound state interpretation of the Luscher 
formula) gives a similar, but somewhat smaller estimate of A/a ~ 0.75 GcV for the lightest 
quark mass. A smaller value of A assumes a slower falloff with box size and thus allows for 
a larger finite size effect on the 16'^ and 20'^ lattices. Thus, in order to obtain a conservative 
upper bound on these effects, we have assumed a value A/a = 0.7 GeV, i.e. slightly smaller 
than the RQM and wave function estimates. To determine the sensitivity of the conclusions 
to the value of A, we also fit the data using A/a — 0.9 GeV, which yields an upper bound on 
the 16^ box about 50% smaller than the A/a — 0.7 GeV fit. In Table H, we give the results 
of fitting the 12^, 16'^, and 20^ LQCD Monte Carlo results to the finite volume formula Eq. 
(p9[). For aEi and Jb, the Table gives the fitted infinite volume result, and an upper bound 
on the finite volume term evaluated on 16'^ and 20^ lattices. 

From Table || it is seen that, in all cases, the estimated finite volume effect on both aEi 
and on /b is smaller than our statistical error by more than a factor of two on the 16'^ lattice 
and by more than an order of magnitude on the 20^ lattice at /3 = 5.9. To determine the size 
of finite volume effects on the quantities Mb, — Mb^ and fB^/fs^, we also need to estimate 
the error on the slope parameters E[ (0) and f'g in Table g. From the K-dependence of the 
finite volume fit parameters, we estimate an approximate upper bound on the finite volume 
error for the slopes to be {AE[)/E[ < .03 and {Afg)/fg < .05 for the 16^ box, and by 
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Table 9: Estimate of finite volume corrections to heavy-light results for (3 = 5.9. 
For each n, results include the fitted infinite volume value for the ground state 
energy and decay constant, as well as estimated upper bounds on the finite volume 
corrections on 16^ and 20'^ boxes. Unbracketed and bracketed numbers result from 
fitting to a Luscher asymptotic form and to a power law (L^^) form, respectively. 



K 


0.154 


0.157 


0.158 


0.159 


Kc 


aEi{oo) 


.719(5) 


.678(6) 


.664(7) 


.645(9) 


.638(9) 




[.719(7)] 


[.679(9)] 


[.664(11)] 


[.638(13)] 


[.636(13)] 


A{aEi){l6) 


± .002 


± .002 


± .003 


± .004 


±.003 




[± .006] 


[± .008] 


[± .008] 


[± .012] 


[±.011] 


A{aEi){20) 


± .0003 


± .0004 


± .0004 


± .0006 


±.0005 




[± .003] 


[± .004] 


[± .005] 


[± .006] 


[±.006] 


/i3(oo) 


.346(11) 


.300(11) 


.283(12) 


.257(14) 


.246(14) 




[.348(16)] 


[.299(17)] 


[.280(19)] 


[.248(23)] 


[.244(21)] 


A/b(16) 


± .004 


± .004 


± .004 


± .006 


±.005 




[± .013] 


[± .014] 


[± .015] 


[± .006] 


[±.005] 


A/b(20) 


± .0006 


± .0007 


± .0007 


± .0011 


±.0008 




[± .007] 


[± .007] 


[± .008] 


[± .011] 


[±.009] 



{AE[)/E[ < .005 and (A/^)//^ < .008 for the 20^ box. Again this is about a factor of two 
below our statistics for 16'^ and entirely negligible for 20"^. 

Recently, it has been argued [Q that, in intermediate ranges of volume where the 
asymptotic behavior predicted by Luscher's volume formula has not yet set in, the volume 
dependence might be expected to exhibit a power law dependence of the form 

const. 



m{L) = 77i(oo) + 



£3 



(42) 



instead of the exponential falloff of Luscher's result. This power law form is also found by 
the authors of Ref. [{48| to fit better to their data on light hadron masses (in full QCD). If 
we assume a similar power-law dependence for the heavy-light data, we obtain extrapolated 
infinite volume results and bounds on finite volume corrections which differ from those 
obtained with Luscher's form. The values in Table which are enclosed in square brackets 
are the results obtained by assuming a power law dependence of the form (H2). Notice 
that the extrapolated infinite volume values change very little compared with the previous 
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analysis. The bound on the finite volume effects at 16'^ are somewhat larger, while those on 
the 20^ lattice are much larger. However, in all cases, the bound on the finite volume effect 
is less than the statistical error. 

To estimate finite size effects for the other (3 values, it is reasonable to assume approx- 
imate scaling. The box sizes for the other /3's have been selected so that they are all of 
about the same physical size as the 16'^ box at /3 = 5.9 (between 1.8 and 2.0 fm). Thus, 
we conclude that finite size effects on all of our data is smaller than our present statistical 
errors. In order to quote a systematic error on our final results (see Section 6.1) for /s^//b„, 
Mb^ — Mb^ and /s, we have assumed that the percentage errors for the other (3 values are 
the same as those obtained at /? = 5.9 on the 16"^ box. 

5.3 Extrapolation to Kc 

To investigate the sensitivity of the chirally extrapolated mass and f-values to the fitting 
range in k, we have done a detailed study of the dependence of the results of the global (in 
k) chi-square fit (^6|) on the k values chosen, for the case /3 = 5.9 on a 16'^ lattice. For this 
run, correlators were studied at k values of 0.154, 0.156, 0.157, 0.158, and 0.159 (with the 
critical Kc =0.15975). The fits were done using a Euclidean time window r< = 3, r> — 10. 
By taking various subsets of k values to perform the chiral fit (cf. discussion at end of 
Section(3)), we can probe the sensitivity of our results to the assumption of linearity of 
mass and f-values in k~^. The central values obtained from the fit (together with the 
associated statistical errors) are displayed in Table O. The range of kappa values used in 
the fit is indicated in the first column using the abbreviated notation n ~ 0.15a; — > x (thus: 
467 indicates that the extrapolation to Kc was made using k values 0.154, 0.156, and 0.157). 



Referring to Table 10, we see that the variation in the extrapolated ground state mass 
and f- value obtained by choosing three sliding windows of adjacent kappa values are in every 
case considerably smaller than the associated statistical errors. For the mass, the central 
values vary by about 50 % of the statistical error, while for the f-value the variation is 20-25 
% of the statistical error. Even totally nonovcrlapping fits (rows 467 and 89) give central 
values lying well within the statistical errors. The statistical errors of course tend to increase 
as we approach Kc ; it is more difficult to detect a systematic trend in the central values 
because the dominant errors are statistical. For the slopes (derivatives with respect to k^^) 
needed for the extraction of Bs properties, the situation is similar. Aside from the 89 fit. 



35 



Table 10: Estimate of systematic effects in chiral extrapolation. Listed are the 
values at Kc and slopes (as a function of k~^ — k~^) of the ground state energy and 
decay constant obtained from a fit to subsets of the five k values, .154, .156, .157, 
.158, and .159. First column indicates the set of k's used by listing the last digit of 
each K included. 



K, range 


Ei{Q) 


Em 


/b(0) 


/s(0) 


X^ldof 


467 


.643(10) 


.319(22) 


.464(29) 


.689(68) 


1.03 


678 


.645(10) 


.303(36) 


.463(29) 


.698(110) 


0.96 


789 


.642(13) 


.332(72) 


.458(36) 


.744(202) 


0.87 


89 


.640(15) 


.385(123) 


.453(41) 


.832(346) 


0.84 


6789 


.643(12) 


.321(52) 


.459(33) 


.725(146) 


0.89 


46789 


.643(11) 


.320(33) 


.460(31) 


.700(95) 


0.92 



which gives a poor determination of the slopes, the central values for all the subsets of k's 
are well within a standard deviation of the full ht to all five kappa values. 

We may conclude from the preceding that, as in the case of finite volume corrections, 
nonlinearities in the chiral extrapolation are not an important source of systematic error in 
our results. In order to arrive at an actual estimate of the chiral extrapolation contribution 
to the total systematic error we have taken the variation in the 3-kappa fits in Table (pXJ) (i.e 
rows labelled 467,678, and 789) which have a reasonable lever arm in 1/k, and fairly small 
statistical errors, as an indication of the extrapolation error to Kc (i-e. in Jb^ and Mb^)- 
As we measure quite close to Ks (at k =0.157), there is effectively no extrapolation error in 
the strange quark quantities. We assume that the chiral extrapolation at [3—b.9 is typical of 
other 13 values. In this way a chiral extrapolation part of the total systematic errors quoted 
in Section 6.1 can be obtained. 

5.4 Scale Errors 

In order to quote physical values of masses and decay constants, one must select a particular 
dimensionful quantity to define the scale. In our discussion, we have taken the values of 
a~^ obtained from the IP-IS charmonium splitting p3{ ai (3 — 5.7,5.9, and 6.1. Our choice 
of a^^ aX (3 — 6.3 is obtained by evolving from (3 = 6.1 via one-loop asymptotic freedom. 
(The same value of a^^ = 3.08 is also obtained from our value for mp.) Other possible 
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choices for the scale-defimng parameter include string tension, rho mass, and /^. Since the 
quoted values of the decay constant /b include a factor of aa , it is particularly important to 
estimate the possible systematic error in our results arising from uncertainty in the overall 



scale at each f3 value. In Fig. 5.4 we have plotted the scale obtained from nip (circles and 
filled circles) , /^r (squares) , and string tension (diamonds and filled diamonds) relative to the 
scales chosen in this paper (Table p. Also included on the plot are points (filled squares) 
obtained from lattice calculations of the dcconfincmcnt temperature Tc [|38| . Since the 
experimental value of Tc is not known, these calculations only give a relative determination 
of the scale at different /3's. (The absolute scale for these points has been chosen to be 
equal to that in Table |l| at /3 = 6.3, which corresponds to a deconfinement temperature of 
kTc = 264 MeV.) The values for mp are from GFll ^^ (filled circles) and from our data 
(Table Q) (circles), while the string tension is a combination of our results at 5.7, 5.9, and 
6.1 (diamonds) and those of Ref. |^ at 6.0, 6.2, and 6.4 (filled diamonds). For the latter 
points, the charmonium scales were estimated from Table H by linear interpolation in In a. 
The values for /^ are taken from Table H. 



The trend exhibited by the data in Fig. 5.4 indicates a significant scale discrepancy 
in the range (3 = 5.7 to 6.0, with the nip scale being about 10-20% higher and the string- 
tension scale about 10-15% lower than charmonium. For f3 > 6.2, the scales appear to 
converge to much better agreement, with deviations of < 5%. This suggests that much of 
the discrepancy at lower (3 is due to finite lattice spacing effects, as opposed to being an effect 
of the quenched approximation. (Discrepancies which do not go away in the scaling limit 
can be ascribed to the neglect of closed quark loops.) The data shown in Fig. |] illustrates 
that, over the entire range of /3, the charmonium scale differs little from a weighted average 
of the other choices. This provides some additional confidence in our choice of scales. To 
estimate the scale error on our heavy-light results, we have used the charmonium scale errors 
quoted in Ref. [Q (which include both statistics and systematics). For /3 — 6.3 we have 
taken a conservative scale error estimate of 5%, based on the spread of values shown in Fig. 



5.4 



Our final results for Mb^ — ^b^ and fs^ {IbJIb^ is dimensionless) quoted below 
are therefore subject to a 5% and 7% error (resp.), assuming that continuum extrapolated 
objects are determined primarily by the larger (3 values where the scale discrepancy is small. 
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The larger error for /b arises from the fact that the quantity computed on the lattice scales 
like a3/2. 

5.5 Extrapolation to the continuum 

By far the largest systematic error in our calculations arises in the extrapolation of the fs 
results to zero lattice spacing. In comparison, the systematic errors incurred from working 
on a finite volume lattice, at finite light quark mass, or even (very probably) the neglect 
of quark loops are negligible. The difficulties here are both intrinsic and practical. On 
the one hand, the detailed form of the lattice spacing dependence of lattice quantities is 
generally rather complicated (involving logarithmic as well as power dependence on the 
lattice spacing) , in contrast to the relatively well understood structure of the chiral or finite 
volume extrapolations. On the other hand, reduction of the lattice spacing by a factor of 
2 requires increasing the lattice volume 16- fold (if we maintain fixed physical space-time 
volume). 



These issues are particularly important in the case of heavy meson decay constants. We 
find that the lattice spacing dependence for fs in the static limit is considerably stronger 



than for fir. This is illustrated in Fig. 5.5, Although our calculations have led to reasonably 
precise results at finite a, quoting a systematic error on the continuum extrapolated result 
clearly requires an investigation of the variation induced by alternative fitting procedures. 



Table 11: Comparison of continuum values for a linear vs. a quadratic fit to the 
a-dependence of physical quantities. Values shown are the fitted a = values of 
/b, Mb^ — Mb^, and /b^/Zb^, and the x^ per degree of freedom for each fit. Errors 
shown are extrapolated statistical errors. 



fit 


fB{MeV) 


xVdof 


Mb, - Mb,, 


xVdof 


fsJfB^ 


xVdof 


Co + cia 
Co + cia^ 


188(23) 
214(13) 


.46/2 
.40/2 


86(12) 

80(7) 


2.4/2 
2.3/2 


1.22(4) 
1.21(2) 


2.4/2 
2.2/2 



We have investigated this variation by fitting the /b results at finite /3 to various func- 
tional forms for the finite a correction. In Table |ll| we show the extrapolated continuum 
result for /b and the x^ per degree of freedom for the alternatives that the finite lattice 
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spacing correction is of order a and of order a^ . The linear and square fits are, explicitly, 

fB{a) = f 88(23) (l + .51(19)a) MeV (43) 

fB{a) = 214(13) [l + [.60(ll)a]^) MeV (44) 

where a is in GeV~^. For Mb^ — Mb^,, the fits are 

Mb,-MbM) = 85.8(11.7) (1 - .24(22)a) MeV (45) 

Mb^-MbSo) = 80.1(6.7) fl - [.47(40)0]^") MeV (46) 



and 



IbJIbSo-) = 1.216(41) (l-.02(6)a) (47) 

IbJIbM = 1.213(23) (l-[.16(31)af) (48) 



These fits are shown with the data in Fig. 5.5. It is clear from this graph that it would 



be very difficult to distinguish between these two possibilities by calculations in the range 
/3 — 5.7 — 6.3, even with improved statistics and additional /3 values. In the absence of more 
precise data at much larger /3 or a complete control of all 0{a) lattice effects, we will take 
the variation of the extrapolated results in Table O as an estimate of the systematic error 
associated with extrapolating to a = 0. In our final result for /b, we quote the continuum 
value obtained from the linear fit with asymmetric errors of +26 and -0 to allow for this 
extrapolation uncertainty. 
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6 Physical Results and Discussion 

6.1 Final results 

In this final section, we collect our results for /_b,//b„, Mb^ — Mb^, and /b^, at each 
lattice spacing calculated and then present the results extrapolated to the continuum. The 
summary of results for each of the four lattice spacings studied are collected in Table R2. A 
comparison with other recently reported results follows immediately In all cases the first 
error quoted is a statistical error obtained by the jackknife procedure described previously. 

Table 12: Final results at fixed (3 (lattice spacing) for Jb, Mbs—M^u, and Jbs/ Jbu in 
the static approximation. The first error is statistical and the second is systematic. 



(3 


/s(MeV) 


Mb, - MBuiMeY) 


Ibs/Jbu 


5.7 


271 ±13 ±20 


66 ± 7 ± 6 


1.181 ±.030 ±.012 


5.9 


241 ±13 ±13 


74 ± 5 ± 6 


1.211 ±.018 ±.014 


6.1 


215 ±21 ±14 


87±9±7 


1.226 ±.027 ±.016 


6.3 


225 ±17 ±14 


68 ± 10 ± 5 


1.172 ±.031 ±.011 



The source of systematic errors varies somewhat with the quantity being computed. For 
./Bs/./Sui the uncertainty in the scale cancels, as the ratio is a dimensionless quantity. Thus, 
for this quantity, the quoted systematic error includes finite volume and chiral extrapolation 
errors. The lattice spacing dependence of this quantity (see Fig(7)) is very small, so we have 
not included a continuum extrapolation error. The result is 



IbJIb^ = 1.216 ± 0.041(stat) ± 0.016(syst) 



(49) 



For the mass splitting Mb^ — Mb^, the systematic errors include finite volume effects, 
the chiral extrapolation (to determine Mb^), and an estimate of the scale error. For the 
continuum extrapolation, we quote the result obtained from the linear fit in a, and take the 
difference between the two fits shown in Table 11 as an estimate of our extrapolation error 
(+0 to —6). Combining the extrapolation error with our other systematic errors, we obtain 
the result 

Mb, - Mb^ = 86 ± 12(stat) t9(syst) MeV. (50) 

Finally, for Jb itself, there is, in addition to the usual finite volume and chiral extrap- 
olation errors, a substantial (not necessarily linear) lattice spacing dependence, as well as 
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the possibility of sizable renormalization corrections to Z at the two and higher loop level. 
These additional systematic errors are quoted separately in our final result: 

/b„ = 188 ± 23(stat) ± 15(syst) If (extrap) ± 14(renorm) MeV. (51) 

6.2 Comparison with other results 

There have been several recently reported studies of the heavy-light meson system in the 
static limit of quenched LQCD. In this Section we will assess the results reported here in 
comparison with these other studies. We will consider four quantities, -Ei, the heavy- light 
ground state energy, the Mb^ ~Mb^ mass splitting, fs the ground-state decay constant, and 
the ratio fs,! Ib^- Even before the present study, there has been some apparent disagree- 
ment among different groups over the size of /s . Some of these discrepancies can be traced 
to different choices of (3 and k, different definitions of the QCD length scale, and different 
evaluations of the perturbative renormalization constant Z . Since our data has provided a 
more accurate determination of the dependence on k and /3, we are able to interpolate our 
data in both variables and make a direct comparison with other groups. 

The focus of this subsection is on the lattice measurement of the ground state energy 
eigenvalue Ei and its matrix element with the unrenormalized local axial current operator 
on the lattice, Jb = vf (c.f. Eq.(l22)). These are the quantities that are extracted directly 
from the lattice heavy-light meson propagators and, for given /3 and k, are independent of 
the choices for length scale and renormalization constant. 

First consider the ground state energy Ei. In Section 4.1 it was shown that the k 
and a dependence of Ei is most easily described by introducing the naive quark mass 



m 



q 



{2fia) ^ — {2fica) ^- The values of Kc for /3 = 5.74 and 6.26 are taken from ^ 



while that at P — 6.0 is taken from ||l^. The Kc values for /? = 5.7,5.9,6.1, and 6.3 are 
from our own data. The values for the scale a~^ are subject to somewhat more uncertainty. 
For 5.7, 5.9, 6.1, and 6.3, the scales we have used are given in Table |l| In order to have 
a reasonably self-consistent set of scales, the remaining values for a^^ are obtained from 
those at 5.7, 5.9, and 6.1 by a simple linear interpolation or extrapolation in In a. Note that 
we are only using the value of a~'^ here to define the naive quark mass iriq. The numerical 
values of Ei and Jb are obtained directly from the lattice propagators and do not depend 
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on choice of scale. Thus, our comparison of data is insensitive to an overall, uniform change 
of scales. The values for the ground state energy Ei (in lattice units), extrapolated to Kc, 



from the various studies are plotted in Fig. 3.2, We conclude from these plots that all the 



data are in reasonable agreement, both on the magnitude of Ei and on its a dependence. 
This is not surprising since, for any reasonably well-chosen smearing function, the value of 
El is obtained fairly unambiguously from the smeared-smeared correlator. Our multistate 
analysis provides a value of Ei with considerably smaller errors than a single-channel analysis 
(note the much smaller error bars for the ACPMAPS results), but the results are completely 
consistent with previous calculations. 

Results from various groups for the AIb^ ~ Mb^ splitting are compared with ours in 



Fig 3.2. Again, the accurate determination of the slope of the ground state energy with 
respect to 1/k is the key to the much smaller error bars shown for the ACPMAPS results. 
However, the results are basically consistent within errors. For this quantity, the lattice 
spacing dependence is very mild. 

Results from various groups for the ratio of decay constants fs^ //s,, are compared with 



ours in Fig 6.2 



In contrast with the situation for Ei or Mg^ — M^^ , there is significant disagreement on 
the value of fg among the various studies, even after correcting for the different values of k 
and p. Results from references [^, 41, IJ, |lj, Q are compared with our data in figure JG^ 



This is a plot of fso-^^ against lattice spacing a, so that we can compare 'raw' lattice 
results without renormalisation factors. The results of [Q (at (3 = 6.3) and |Q are in good 
agreement with ours, while [|^ ^ Q report substantially larger values. 

We believe that the discrepancies can be traced primarily to the incomplete isolation 
of the ground state with the sources chosen. Observe in particular that the result in 149], 
where a variational method is used to isolate the ground state, is in good agreement with 
us. Sensitivity of the apparent value of Jb to the form of the source has been emphasized 
inisi. 
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6.3 Future Studies 

We have presented results for the decay constant fg and for masses of low-lying heavy-light 
states in the static approximation. The analysis procedure introduces several improvements 
over previous smearing methods. First, the success of the RQM in reproducing the measured 
lattice wave functions is exploited by using the RQM to construct not only an accurate 
ground state wave function, but also a set of orthonormal excited state smearing functions. 
Second, we make full use of the information contained in the matrix of smeared-smeared and 
smeared-local correlators, including both ground state and excited state smearing functions 
at each end. Our method provides much greater control over systematic errors from higher 
state contamination, because of the fact that the source smearing functions arc tuned directly 
to the lattice wave functions, without regard to the behavior of the effective mass plots. The 
appearance of long plateaus in the SS and SL plots at the same value of effective mass is thus 
strong evidence that the systematic error from higher states has been largely eliminated. 
We will report the application of our methods for heavy-light mesons to the spectrum of 
radial and orbital excitations for heavy-light systems and the B parameter in forthcoming 
papers. 

We expect to be able to improve the accuracy of the present results for the Mb^ — 
Mb^ and fs^/ Jb^ by using larger ensembles. There are two other sources of systematic 
uncertainty in our results for Jb- Use of the Wilson action for the light quarks implies 
lattice spacing corrections in 0(a) and the large one loop rcnormalization for the axial 
current suggests that the two loop correction may be sizable. Study of these effects will be 
required to substantially improve the error on Jb- 
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A Light Hadron Results 

In this Appendix we report the resuhs for the light quark sector which we have used to set the 
parameters for the heavy-hght meson properties. First consider the light-light pseudoscalar 
meson mass as a function of the light quark hopping parameters. This serves to determine 
the critical hopping parameter at each value of [3 studied. Using the physical K mass, we 
also establish the correct hopping parameter for the strange quark. (The small effect from 
the nonzero up or down quark mass has also been included in our results.) In addition to 
determining the light quark hopping parameters from the pseudoscalar meson masses, the 
meson propagators have also been analyzed to obtain the values of /^r and nip. All of the 
results discussed here for the light hadron parameters have been extracted from local qq 
operators. An analysis of light meson parameters using smeared sources in Coulomb gauge 
is in progress and will be reported elsewhere. 

Let the hopping parameters for the two valence quarks in the light meson be denoted 
by Ki and K2, and denote the pseudoscalar meson mass by mp. All of our data is consistent 
with a linear dependence of m?p on k^ and k,2 ■ The results can be expressed in terms of 
the parameters C and k^ of the linear fit 

(mpa)2 =C(kj;1 + k^1-2k-i) (52) 

The results for the four [3 values studied are given in Table |3[ Using the scales in Table |l|, 
the naive quark masses, defined by 

1 1 

(53) 



2K,a 2Kca 

are also listed in Table O. (Here, up and down quarks are taken to be degenerate in mass.) 
A plot of TOp vs. quark mass (defined in Eq. (p3|)) is shown in Fig. ^. For [3 — 5.7, 6.1, and 
6.3 the values of mpa were obtained using equal quark masses, ki — K2. For /3 = 5.9 several 
combinations of unequal quark mass values (open circles in Fig.H) were used to check the 
validity of Eq.(|2|). For f3 = 5.9, in addition to the main run on a 16^ x 32 lattice, in order 
to investigate finite volume effects, we have also carried out Monte Carlo runs on 12'^ x 24 
and 20^ x 40 lattices. The results labeled 5.9(oo) are the infinite volume values obtained 
by fitting all three box sizes to the functional form derived by Luscher, as discussed in Sec. 
5.2. 
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Table 13: Light quark hopping parameters obtained from Hght pseudoscalar meson 
mass measurements. Values are shown for each run, as well as the infinite volume 
extrapolated value at /? = 5.9. C is the slope of the pseudoscalar (mass)^ as a 
function of k^^. The last two columns give the strange and up quark masses x2a. 
Errors are statistical. 



Run 


/3 


Kc 


C 


^s '^c 


'^u '*c 


h 


5.7 


.16914(10) 


.703(10) 


.2533(36) 


.01024(15) 


e 


5.9 


.15972(14) 


.615(14) 


.1209(27) 


.00489(11) 


c 


5.9 


.15975(6) 


.609(9) 


.1221(18) 


.00494(7) 


f 


5.9 


.15981(4) 


.591(8) 


.1258(17) 


.00509(7) 




5.9(oo) 


.15980(4) 


.597(7) 


.1245(14) 


.00504(6) 


d 


6.1 


.15496(3) 


.480(11) 


.0832(19) 


.00336(8) 


g 


6.3 


.15178(4) 


.395(10) 


.0613(15) 


.00248(6) 



In Table |l^, the values for C and Hc are independent of the scale chosen for a^^. The 
remaining columns are computed using the scales in Table y. The errors quoted in Table |l3| 
are statistical only, and do not include the uncertainty associated with the choice of a^^. 
The uncertainty arising from the scale determination, along with other sources of error, will 
be discussed in Section 6. 

Although they are not used in the body of the paper, we also briefly discuss our results 
for the rho mass and the pion decay constant. We again emphasize that all of the results 
discussed here are obtained from light quark propagators with (5-function sources. The 
results are thus subject to possible systematic errors from higher state contamination. This 
is not a problem in measuring the pion and kaon masses discussed above, but it becomes 
more of a difficulty for mp and /^ measurements. In order to determine /^, one must 
calculate both the propagator with a pseudoscalar (9759) source at each end (PP), and the 
propagator with a pseudoscalar source at one end and an axial vector (^75709) source at the 
other end (PA). (Note: One may also use the (AA) propagator with an axial vector source at 
both ends. However, since the vacuum to one pion matrix element of the axial vector source 
contains an explicit factor of 777,^, the pion pole residue in the (AA) propagator vanishes 
more rapidly than that of the (PA) propagator in the chiral limit, making it more difficult 
to measure accurately.) The (PP) and (PA) propagators are fit to a single exponential in 
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Table 14: Results for light-light mesons. Meson propagators were fit to a sin- 
gle exponential over time window AT. Result for /^r includes a tadpole improved 
renormalization factor, computed in Ref.[20]. 



p 


AT 


K 


TTiT^a 


nipa 


ha 


5.7 


6-10 


.161 


.649(3) 


.787(2) 


.145(3) 






.165 


.456(5) 


.675(5) 


.133(4) 






.1667 


.351(8) 


.629(15) 


.115(6) 






.168 


.237(13) 


.586(44) 


.102(11) 






Kc 





.566(11) 


.108(14) 


5.9 


8-12 


.154 


.527(2) 


.619(2) 


.104(2) 






.156 


.426(3) 


.546(3) 


.093(2) 






.157 


.360(2) 


.513(3) 


.086(2) 






.158 


.288(3) 


.479(5) 


.078(2) 






.159 


.195(4) 


.444(12) 


.072(4) 






Kc 





.418(4) 


.067(2) 


6.1 


10-16 


.151 


.409(3) 


.482(3) 


.081(2) 






.153 


.276(4) 


.401(4) 


.069(2) 






.154 


.196(3) 


.361(8) 


.061(2) 






.1545 


.137(5) 


.341(19) 


.056(4) 






Kc 





.324(8) 


.055(2) 


6.3 


12-18 


.1500 


.249(3) 


.331(3) 


.054(2) 






.1510 


.160(6) 


.283(5) 


.045(4) 






.1513 


.126(6) 


.269(7) 


.041(4) 






.1515 


.099(5) 


.260(11) 


.038(11) 






Kc 





.246(8) 


.037(3) 
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Table 15: Light-light meson results at /3 = 5.9 for lattice sizes 12^ x 24(e), 16^ x 32(c), 
and 20^ x 40(/). Result for /„ includes the tadpole improved renormalization factor, 
computed in Ref. [20]. 



p 


AT 


K 


m^a 


JTlpfl 


A a 


5.9(e) 


8-12 


.154 


.535(4) 


.618(4) 


.105(3) 






.157 


.364(4) 


.505(7) 


.083(4) 






.158 


.306(9) 


.475(15) 


.080(4) 






.159 


.198(33) 


.375(77) 


.075(17) 






Kc 





.403(14) 


.067(5) 


5.9(c) 


8-12 


.154 


.534(3) 


.622(2) 


.106(3) 






.157 


.364(4) 


.508(5) 


.087(3) 






.158 


.292(4) 


.469(8) 


.079(4) 






.159 


.188(8) 


.440(27) 


.067(8) 






He 





.407(8) 


.067(4) 


5.9(f) 


8-12 


.154 


.525(2) 


.617(2) 


.103(2) 






.157 


.358(3) 


.515(4) 


.085(2) 






.158 


.288(3) 


.482(5) 


.078(3) 






.159 


.196(4) 


.442(12) 


.073(4) 






Kc 





.423(5) 


.068(2) 
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the time ranges shown. Table W summarizes our results for light hadron parameters. Again, 
the errors quoted are purely statistical. Some systematic errors are expected particularly 
for /tt, for which a stable mass plateau in the (PA) propagator with the same mass as 
the (PP) propagator was not generally achieved. An attempt was made to compensate for 
this by fitting the (PA) propagator with the mass fixed to be equal to that of the (PP) 
propagator. It is clear from these results that a much better determination of /^r from our 
data will be possible when smeared operators are employed. The rho propagators exhibit 
reasonable plateaus in the time intervals shown in Table |A[ but a study of the variation of 
the results with different AT windows indicates that a systematic error of from 1 to 2 times 
the statistical errors cannot be ruled out. The values of f^a given in Table ^ include the 
perturbative renormalization constants computed in Ref. pQ]. For (3 = 5.9, the results in 
Table ^ are the infinite volume extrapolated results from the three Monte Carlo runs on 12"^ 
(e), 16'^ (c), and 20"^ (f) boxes (except for k — .156 which was done only on the 16"^ box). 
The results on each size box are listed separately in Table H. In the main analysis of this 
paper, the only light meson parameters we will use are the hopping parameters determined 
from the pseudoscalar masses. Since the (PP) propagator always exhibits a stable mass 
plateau, these parameters are well-determined and should be relatively free of systematic 
error from higher states. 
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B Multistate Extraction of Meson Wavefunctions 

A particularly graphic illustration of the power of the smearing technique in reducing 
the contamination of higher states is obtained by examining the time development of the 
Coulomb gauge Bethe-Salpeter wavefunction of a static-light meson beginning with either 
(a) the RQM smeared source defined in equation (|l7|), or (b) a cube smeared source |4l| . 
In the first case, we extract 

*B5(i?,r) =< 0|g(i?,r)Q(0,T)|$('^\0 > (54) 

for the ground state (a—l) at small Euclidean times. We have done the comparison for the 
case /3 =5.9, k =0.159, on a 16'^ lattice. In Fig.H the evolution of the wavefunction using 
a source smeared with the ground state wavefunction of the RQM (with constituent mass 
/i=0.12) is shown for Euclidean times T=l,2 and 4. It is apparent that the wavefunction 
has reached its asymptotic value to very good accuracy already at time slice 2, with little 
further change at T—i (in fact, the overlap of the wavefunction at T—2 with that at T—A 
is 0.9986!). In Fig.^ the corresponding evolution (again for Euclidean times T—1,2 and 4) 
is shown for a source smeared over a cube of width 7 lattice spacings. Here the convergence 
is much slower, with the pointwise convergence near the origin particularly tardy. 

Although the use of a single smearing function obtained from the RQM is adequate 
to the task of extracting the ground-state Bethe-Salpeter wavefunction, even the improved 
smearing given by the RQM is not sufficient if we wish to do the same for the higher excited 
states in a given channel. In Section 4 we showed how to define optimized smeared states 
I ^'^ > in which the admixture of all but one (the A'th) of the first M meson states in 
a given channel is tuned to zero. For example, taking M=3 at /3=5.9, k=0.159 on a 16'^ 
lattice, one finds that the choice of smearing function (cf. Sect 3.2,Eq (p3)) 

I $(2) >= 0.03 I $(^) > -fO.82 I $(2) > ^0.57 | $(^) > (55) 

produces an optimized first excited state in the sense that admixtures of the ground and 
second excited state are tuned out (based on a fit of the form (P2|)) with r<=2,T'>=7). The 
use of such an optimized smearing is crucial if we wish to extract the correct Bethe-Salpeter 
wavefunction of the first radial excited state near the origin. Any sizable admixture of the 
exact meson ground-state will otherwise dominate the small r region of the wavefunction at 
large Euclidean time, before the higher (2nd, 3rd, etc) states have decayed away. In Figure 
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H wc show this phenomenon with a Bcthe-Salpeter wavefunction ^B5(i?,r) defined as 

^Bs{R,T) =< I q{R,T)Q{0,T) I $(2)^0 > (56) 

The wavefunction (renormaUzed to unit norm) is plotted for EucHdean times T—1,2,A, and 
6. There is a steady upward drift of the wavefunction at the origin as T increase to 4, but 
by time shce 6 the influence of the ground state is clearly apparent as the latter begins 
to dominate the evolved meson state. It would clearly be very difficult to draw any firm 
conclusions about the behavior of the excited state wavefunction close to the origin from 
these measurements. 

On the other hand, using the optimized smearing found above (pq), and computing 

^TsiR^T) =< I q{R,T)Q{0,T) I ¥^\0> (57) 

one finds (Fig.H shows the optimized wavefunction for times T— 1,2,3) a rapid convergence 
to an asymptotic shape by the third time slice, giving a value for the wavefunction at the 
origin ~0.44, as compared to a maximal value ~0.36 obtained at T=4 from Eq ( pq ) before 
convergence is lost. Eventually, of course, the ground state will dominate in this case also, 
but by using the optimized state, we correctly extract the exact Bethe-Salpeter wavefunction 
for the excited state before the ground state contamination has a chance to become sizable. 
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Figure Captions 
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Figure 1: Static potential calculated in Coulomb gauge at /3 = 5.7,5.9,6.1 and, 6.3 on 
lattices of size 12'' x 24, 16^ x 32, 24^^ x 48, and 32^ x 48, respectively. Errors shown 
are statistical. 

Figure 2: Smeared-smeared (S-S) and smeared local (S-L) effective masses for /? — 5.7, 12^ x 
24, and n = .161, .165, .1667, and.lGS. Solid line is the ground state energy extracted 
from a 2-state fit. Smearing functions are optimized combinations of RQM wavefunc- 
tions, as described in the text. 

Figure 3: Smeared-smeared (S-S) and smeared local (S-L) effective masses for /3 — 5.9, 12"^ x 
24, and n = .154, .157, .158, anrf. 159. Solid line is the ground state energy extracted 
from a 2-statc fit. Smearing functions are optimized combinations of RQM wavefunc- 
tions, as described in the text. 

Figure 4: Smeared-smeared (S-S) and smeared local (S-L) effective masses for f3 — 5.9, 16"^ x 
32, and k = .154, .157, .158, and. 159. Solid line is the ground state energy extracted 
from a 2-statc fit. Smearing functions are optimized combinations of RQM wavefunc- 
tions, as described in the text. 

Figure 5: Smeared-smeared (S-S) and smeared local (S-L) effective masses for (3 — 5.9, 20^ x 
40, and k = .154, .157, .158, and. 159. Solid line is the ground state energy extracted 
from a 2-state fit. Smearing functions are optimized combinations of RQM wavefunc- 
tions, as described in the text. 

Figure 6: Smeared-smeared (S-S) and smeared local (S-L) effective masses for /3 — 6.1, 24"^ x 
48, and k = .151, .153, .154, an(i.l545. Solid line is the ground state energy extracted 
from a 2-state fit. Smearing functions are optimized combinations of RQM wavefunc- 
tions, as described in the text. 

Figure 7: Smeared-smeared (S-S) and smeared local (S-L) effective masses for /3 — 6.3, 32"^ x 
48, and k = .1500, .1500.1513, and. 1515. Solid line is the ground state energy ex- 
tracted from a 2-state fit. Smearing functions are optimized combinations of RQM 
wavefunctions, as described in the text. 

Figure 8: Heavy-light ground-state energy aEi vs. bare quark mass for (3 = 5.7,5.9,6.1, 
and 6.3 Each data point is the ground state energy extracted from a two-state fit. 
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Error bars are statistical only. Solid lines are obtained from a simultaneous two-state 
fit to all kappa values for a given (3 as described in the text. 

Figure 9: Jb as a function of bare quark mass for the four runs /? = 5.7,5.9,6.1, and 
6.3. Data points are decay constants extracted from a two-state fit. Error bars are 
statistical only. Solid lines are obtained from a simultaneous fit to all kappa values 
for a given beta as described in the text. 

Figure 10: Heavy-light ground-state energy at k — Kc vs. lattice spacing for j3 — 
5.7,5.9, .6.1, and6.3. Solid line represents a minumum-x^ linear fit to the four data 
points. 

Figure 11: Figure 11: Mb^~Mb^ vs. lattice spacing. Solid line is the best linear fit. Dashed 
line is a quadratic fit (a^) used to estimate systematic error in a ^ extrapolation. 
(See Section 5.5). 

Figure 12: /bs/Ibu vs. lattice spacing. Solid line represents the best linear (a) fit. Dashed 
line is a quadratic fit (a^) used to estimate systematic error in a ^ extrapolation. 
(See Section 5.5). 

Figure 13: Time-window dependence of effective mass for 1-state fits compared with 2-state 
fits for P ~ 5.7, K = .161. The four sets of points represent the different values of 
RQM mass parameter /i used to construct smearing functions. One-state fits were 
obtained from time-windows (reading from right to left on the graph) 1-6, 2-7, 3-8, 
4-9. Points are plotted at exp[— (A_E)t„ii„], where AE = E2 — E1 is the splitting of the 
first excited (2S) state from the ground-state and tmin is the smallest time included 
in the fit. The 2-state fits were obtained from the window 2-8 but are plotted here at 
< « 00 to illustrate the convergence of the 1-state effective mass to the essentially /i 
independent 2-state value. 

Figure 14: Time-window dependence of /^ for 1-state fit and comparison with 2-state fit. 
(See item of Fig. 13.) Also plotted for comparison are the values obtained by ex- 
trapolating the 1-state fits to t = 00, using the measured value of the energy-splitting 
E2~Ei^ 0.32. 

Figure 15: Scales obtained from rup (circles and filled circles), /^r (squares), string tension 
(diamonds and filled diamonds), and deconfincment temperature Tc (filled squares) 
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relative to the scales in Table 1. Our data is denoted by open symbols. 

Figure 16: /s at k = Kc as a function of lattice spacing. The scales and renornialization 
constants used are given in Table 1. The solid line is the best linear fit. The dashed 
line is a quadratic (a^) fit used to estimate systematic error in a ^ extrapolation. 
(See Section 5.5). 

Figure 17: Ground state energy, in lattice units, at n — Hc, versus lattice spacing. Points 
at the same a have been slightly displaced for readability. 

Figure 18: Mass splitting, Mb^ — Mb^, in physical units, versus lattice spacing. Points at 
the same a have been slightly displaced for readability. 

Figure 19: Comparison of Jb^/Ibu vs. lattice spacing for present results and other recent 
works. Scale set from all ACPMAPS values, using empirical fit, a^^ linear in j3. 
Points at the same a have been slightly displaced for readability. 

Figure 20: Comparison of Jb vs. lattice spacing for present results and other recent works. 
Points at the same a have been slightly displaced for readability. 

Figure 21: m^ vs. naive quark mass for /3 = 5.7,5.9,6.1, and 6.3. The vertical scales for 
5.9, 6.1, and 6.3 are offset by multiples of 0.2 for display purposes. Points labeled 
5.9e and 5.9u represent mesons with equal and unequal quark masses, respectively. 
Values with unequal quark mass are plotted at the average mass {mi + ?7i2)/2. 

Figure 22: Time evolution of the ground state heavy-light wavefunction for (3 — 5.9, k = 
.159, using a source smeared with the ground state wavefunction of the RQM with 
constituent mass parameter /i = 0.12. 

Figure 23: Time evolution of the ground state heavy- light wavefunction for /? ~ 5.9, k = 
.159, using a source smeared over a cube of width 7 lattice spacings. 

Figure 24: Time evolution of the first excited state (2S) state heavy-light wavefunction for 
(3 = 5.9, K — .159, using a source smeared over the 2S RQM wavefunction. 

Figure 25: Time evolution of the first excited (2S) state wavefunction using a 3-state fit and 
selecting the optimized combination of RQM smearing functions tuned to be orthogo- 
nal to the ground state and second excited state, a Bethe-Salpeter wavefunction and 
optimized smearing. 
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